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ABSTRACT 

The  objective  of  this  work  was  to  determine  the  expected 
deflections  prior  to  and  at  the  rupture  load  for  a  panel  of  hull  plating 
of  an  ice-strengthened  ship.   An  idealized  mathematical  model  of  a  pin- 
nacle of  ice  being  forced  against  the  plating  was  used  as  the  assumed 
loading  mechanism.   The  panel  of  hull  plating  was  modeled  as  a  clamped 
rectangular  plate.   The  idealized  ice-loads  were  defined  as  a  four  sided 
pyramid  of  ice  which  when  forced  against  the  center  of  a  rectangular 
plate  crushed  at  a  given  pressure.   Thus  the  size  of  the  load  bearing 
surface  continually  increased  as  the  applied  force  was  increased, 

A  method  of  analysis  vjas   established  which  provides  a  step-by- 
step  solution  of  the  large  deflections  of  a  plate  and  a  means  of  deter- 
mining the  minimum  expected  failure  load.   The  analysis  is  begun  at  the 
point  V7here  the  load  on  the  plate  from  the  pyramid  of  ice  is  large  enough 
to  cause  plastic  collapse  by  a  three  hinge  mechanism,   A  solution  method 
for  the  deflection  of  a  uniformly  loaded  membrane  forms  the  basis  for  the 
calculations  of  the  plate  deflections.   Then  the  energy  absorbed  by  the 
plate  at  the  point  of  expected  failure  is  calculated  from  an  equation 
that  depends  on  the  total  strain  at  the  point  of  failure. 

The  possibility  of  failure  of  the  plate  at  any  given  load  and 
deflection  may  be  determined  from  the  failure  criterion.   This  criterion 
was  established  by  forming  a  relation  between  a  parameter  representing 
the  geometric  properties  and  a  param.eter  representing  the  material  pro- 
perties of  the  plate.   The  relation  was  formed  from  published  data  on 
plates  that  had  been  tested  to  failure.   The  calculated  failure  parameters 
for  these  plates  were  used  to  form  a  linear  equation  V7hich  correlates 
all  the  experimental  failure  loads.   The  criterion  was  then  checked  by 
comparing  the  results  given  in  other  experiments  on  the  large  deflections 
of  plates. 

The  method  of  analysis  was  then  applied  to  plates  similar  to 
the  hull  plating  of  ice-strengthened  ships.   The  results  showed  that  the 
method  of  analysis  can  be  used  with  a  reasonable  degree  of  accuracy  to 
calculate  the  large  deflections  of  a  plate  subject  to  a  lateral  load  and 
to  predict  the  load  and  deflection  at  which  failure  may  be  expected. 
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NOMENCLATURE 

F,Fj^,Fq  --(lbs)    The  force  behind  the  pyramid  of  ice 

P^  --(psi)    Ice-crushing  pressure 

q  --(psi)    Any  uniformly  distributed  load  over  a 

previously  defined  bearing  area 

M  --(in-lb)    The  bending  moment  about  the  longitudinal 

axis 

--(in-lb)    The  bending  moment  at  the  plastic  collapse 
load, 

--(psi)    The  yield  stress  for  a  material 

--(psi)    The  ultimate  stress  of  a  material 

--(psi)    Equivalent  yield  stress  of  material  in  an 
elemental  strip  of  plate 

--(psi)    In-plane  membrane  stress 

--(psi)    Total  in-plane  stress  due  to  the  membrane 
force  and  bending  moment 

--(lb/in)    In-plane  membrane  force  per  unit  length  of 
plate  edge 

-- (in-lb)    Energy  per  unit  volume  absorbed  by  the  plate 
in3 

R  -- (in-lb)    Total  restraint  energy  per  unit  volume  of  a 

j^^3  material  in  simple  tension 

6y  --(in/in)    The  ultimate  strain  of  a  material  in  simple 

tension  corresponding  to  (J 

6  J.  --(in/in)    The  total  strain  at  a  point  due  to  bending 

and  stretching 

6  --(in/in)    The  strain  in  the  x-direction  at  a  given 

point  in  the  plate 

6y  --(in/in)    The  strain  in  the  y-direction  at  a  given 

point  in  the  plate 

6,  --(in/in)    The  in-plane  strain  of  the  outer  fibers  of 

the  plate 


Mo 

<^o 

O^u 

^p 

'^s 

0-t 

s.s„,s. 

U 

Z,Z.,Z   ,  --(in)    Partial  center  deflection  as  determined 

by  membrane  formula 

w  --(in)    Total  plate  deflection  at  any  point  (x,y) 

WjW.jW   ,  --(in)    Total  center  deflection  for  a  given  load 

T  (radians)   Angle  between  horizontal  and  midplane  of 

plate  at  a  given  deflection 

6  (radians)    Angle  of  rotation  due  to  bending  of  a  plane 

cross-section  through  the  thickness  of  the 
plate 

o   =  w/A  a  non-dimensional  failure  parameter  for 

plate  geometry 

0  =  U/R  a  non-dimensional  failure  parameter  for 

plate  material  properties 

k  Coefficient  within  the  membrane  formula  that 

is  dependent  on  x  and  y 

k.  Coefficient  within  the  derivative  of  Z   that 

is  dependent  on  x  and  y 
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INTRODUCTION 


Discussion  of  Proble 


m 


The  solutions  to  two  interrelated  problems  are  sought  in 
this  analysis.   The  failure  of  clamped  rectangular  plates  may  be 
separated  from  ice-loading  of  plates,  but  the  tx-70  studies  are  brought 
together  here  for  the  consideration  of  a  special  problem  in  ship's 
structures.   The  full  analysis  of  hull  plating  subjected  to  the  various 
loading  conditions  associated  with  navigation  through  areas  where  ice 
is  present  brings  the  two  topics  together.   The  first  problem  is  to 
determine  the  relation  between  the  hull  plating  and  the  ice  that  may 
be  forced  against  the  plating.   The  second  problem  is  encountered  as 
the  force  driving  the  ice  against  the  ship  is  increased  and  the  hull 
plating  deforms  until  plate  failure  or  rupture  occurs.   The  total  pro- 
cess may  be  considered  continuous,  starting  with  the  initial  contact 
with  the  ice,  subsequent  ice  crushing  and  plate  deflection,  and  finally, 
extreme  plastic  deformation  of  the  plate  and  failure. 

The  neod  for  this  analysis  is  generated  more  from  the  design 
of  merchant  ships  which  are  to  operate  in  ice  filled  waters  than  from 
icebreaker  design.   Icebreakers  are  generally  very  conservatively 
designed  because  they  must  be  able  to  withstand  ice-loadings  of  almost 
any  magnitude  during  their  lifetime.   On  the  other  hand,  merchant  ships, 
even  those  with  icebreaker  bows,  have  to  balance  their  design  between 
profit  and  repair  costs.   Their  carrying  capacities  are  greatly  affected 
by  the  addition  of  thicker  steel  plates  along  the  ice-belt.   Furthermore, 
if  there  is  a  ch.iuce  that  the  ship  could  lose  part  of  its  cargo  through 
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a  ruptured  plate,  the  economics  of  the  situation  increases.   Even  more 
important  though,  is  the  rising  concern  over  the  possibility  of  oil 
pollution  of  the  arctic  seas.   In  order  to  prevent  oil  spills,  it  has 
been  proposed  that  a  tanker  be  designed  with  outer  buffer  tanks  filled 
■with  water.   The  outer  hull  would  be  of  thinner  plate  than  the  usual 
single  hull  ship.   The  separation  of  the  two  hulls  would  be  based  on  the 
maximum  deflection  that  could  be  expected  from  the  outer  plating  when  it 
failed  under  ice- loads. 

Idealization  of  Ice-Loads 

The  ice-loads  to  be  considered  will  be  idealized  with  regards 
to  their  geometry  and  point  of  application.   The  present  design  practices 
are  based  on  static  loading  of  the  hull  plates.   The  model  loads  used 
here  will  be  considered  continuous  but  applied  in  incremental  static  steps 
Icebreaker  and  ice-belt  plate  scantlings  are  usually  obtained  by  using  a 
uniformly  distributed  pressure  equal  to  the  ice  crushing  pressure.   Since 
the  actual  ice  crushing  pressure  is  not  constant  and  dependent  upon 
geographical  location,  age,  salinity  and  many  other  factors,  it  was 
desired  to  try  to  cover  a  wide  range  of  possible  loading  mechanisms. 
The  driving  force  behind  the  plate/ice  interaction  must  also  be  considered 
It  was  for  these  reasons  that  a  pyramid  shaped  pinnacle  of  ice  was  chosen 
as  the  idealized  ice-loading  mechanisi 


;m. 


Type  of  Analysis 

Plate  failure  analysis  lies  between  the  theory  of  plasticity 
and  a  handful  of  experimental  tests  that  have  carried  applied  loads  up 
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to  plate  rupture.   Theoretically,  it  would  be  desirable  to  use  the  theory 
of  plasticity  to  analyze  the  plate  deflection  and  failure.   This  was 
considered,  but  after  reading  other  works  in  the  field  and  considering 
the  accuracy  of  the  possible  results,  it  was  decided  to  use  a  more 
empirical  approach.   The  final  analysis  will  be  a  combination  of  plastic 
analysis  and  engineering  judgment.   This  combination  of  approaches  to 
the  problem  allows  a  study  of  the  plate/ice  interaction  within  a  limited 
time  period  with  sufficient  accuracy  of  the  results. 

Load/Deflection  Relations 

Most  of  the  experimental  work  on  rectangular  plates  with  large 
deformation  has  been  done  under  uniformly  distributed  loads.   To  obtain 
any  data  on  concentrated  loads,  it  was  necessary  to  go  back  to  a  report 
of  experiments  conducted  in  1937.   The  plate  deflections  and  stresses 
were  kept  in  the  elastic  range  in  these  tests.   The  results  did  provide 
a  good  starting  point  for  the  shapes  of  the  plates  during  plastic  defor- 
mation under  concentrated  and  partial  loads.   Some  of  the  results  and 
observations  on  circular  plates  under  concentrated  loads  were  borrowed 
for  the  analysis. 

It  was  necessary  to  consider  the  concentrated  or  partial  loads 
in  order  to  present  the  complete  model  of  plate/ice  interaction.   The 
crushing  process  of  ice  in  contact  with  the  ship's  hull  will  begin  as 
a  concentrated  load.   The  ultimate  strength  of  sea  ice  in  compression 
is  usually  low  enough  in  relation  to  the  design  pressures  of  hull  plating 
such  that  the  ice-loads  are  uniformly  distributed  before  permanent  defor- 
mation of  the  plate  begins.   However,  since  the  ice  crushing  process  will 
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require  some  partial  loading,  large  plate  deformations  may  occur  under 
smaller,  but  more  concentrated  loads.   Therefore,  partial  loading  of  the 
plate  will  be  analyzed  in  order  to  consider  all  possible  extremes  for 
the  failure  analysis  of  the  plate. 

Failure  Analysis 

Experimental  results  of  rectangular  plate  failures  are  usually 
an  added  column  of  data  to  the  experimenter's  primary  reason  for  testing. 
The  beginning  of  investigations  into  large  plastic  deformation  and  failure 
is  often  attributed  to  Bach.   His  experiments  conducted  in  1908  formed 
the  basis  for  Hovgaard's  20%  permanent  deflection  rule.    These  tests, 
like  most  of  the  data  obtained  for  rectangular  plate  failure,  were  for 

uniformly  distributed  loads.   Often  the  conEnent  that  the  plate  failed  at 

(T4) 
a  given  pressure  is  the  only  amount  of  data  given.   Greenspon,       in 

(E2  E3) 
evaluating  experiments  by  J.  W.  Day,  ^   '     attempted  to  correlate  the 

failure  data  so  that  plate  rupture  under  a  uniformly  distributed  load 

could  be  predicted.   His  recommendation  was  to  establish  a  limiting 

maximum  deflection  to  span  ratio.   Using  a  statistical  approach,  it  could 

be  stated  that  a  uniform  lateral  load  giving  a  deflection  to  span  ratio 

greater  than  0.10  could  be  considered  approaching  the  plate  failure 

load.   The  only  other  author  who  discussed  plate  failure  to  any  extent 

(El) 
was  Loeser. 

Most  of  the  experiments  as  reported  by  the  authors  listed  in 

the  bibliography  under  Plate  Deflection-Experimental  were  helpful  in 

establishing  a  failure  criterion.   It  was  decided  early  in  the  analysis 

that  an  energy  relation  would  be  used.   This  enabled  correlation  of  the 


1    See  reference  T5 ,  pp.  120-125. 
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failure  data  that  was  available.   Not  all  of  the  data  V7as  useful,  since 
all  material  properties  must  be  known  in  order  to  properly  correlate  the 
experiments.   Since  Day  and  Loeser  included  all  m.aterial  properties  in 
their  reports,  the  results  of  their  experiments  were  relied  on  most 
heavily. 

Outline  of  Procedure 

The  procedure  of  investigation  that  follows  will  continue  in 
a  similar  format  as  used  in  the  introduction.   First,  the  model  of  the 
ice-loads  will  be  formed.   Then,  the  plate  deflection/load  relations 
will  be  established.   At  each  increment  of  loading  or  deformation  of  the 
plate,  the  amount  of  energy  absorbed  by  the  plate  material  should  be 
calculated.   The  relations  for  this  segment  of  the  analysis  will  depend 
on  the  load/deflection  calculations  and  will  serve  to  relate  them  to  the 
failure  criterion.   The  final  step  in  the  analysis  will  be  the  estab- 
lishment of  the  failure  criterion.   This  will  complete  the  model  of  the 
ice/plate  interaction. 

Ordering  of  the  Bibliography 

It  has  been  previously  mentioned  that  this  analysis  falls 
between  the  basic  theory  of  plasticity  and  empirical  relations.   Many 
of  the  assumptions  and  choices  of  methods  are  based  on  the  experiments 
and  observations  of  other  authors.   Each  phase  of  the  analysis  can  be 
performed  by  other  techniques.   However,  within  the  time  limits  of  the 
investigation,  all  cases  were  carefully  studied  so  as  to  make  use  of  the 
best  method  available.   For  these  reasons,  it  was  felt  that  the  list  of 
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references  should  reflect  the  degree  of  influence  of  each  work  upon  the 
final  analysis.  Therefore,  the  references  within  each  subject  area  are 
listed  in  approximate  descending  order  of  influence  or  use. 
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PROCEDURE  FOR  CALCULATING  ICE /PLATE  INTERACTION 
Initial  Assumptions  for  Ice-Loads 

A  four  sided  pyramid  of  ice  as  sliOTvn  in  Figure  1  will  be  the 
loading  mechanism.   The  manner  of  loading  will  be  such  that  it  may  be 
considered  to  be  static  at  all  times. 

The  Included  angle  of  the  pyramid  frustum  will  not  effect  the 
load  bearing  area  of  the  ice  on  the  plate.   Also  it  is  small  enough  to 
cause  only  the  crushing  portion  of  the  ice  to  make  contact.   That  is, 
the  sloped  sides  of  the  cone  will  not  come  in  contact  with  the  plate  at 
any  t  ime . 

The  load  bearing  surface  will  conform  to  the  curved  shape  of 
the  plate  as  if  the  plate  in  the  region  of  the  load  is  subjected  to  a 
uniformly  distributed  load.   Thus,  load  edge  stress  concentrations  will 
be  neglected.   Also,  at  no  tliae  will  the  area  of  load  beariiig  surface  be 
less  than  twice  the  thickness.   Therefore,  shear  forces  will  be  neglected, 

The  load  bearing  area  will  be  a  square  of  side  equal  to  C. 
If  C  equals  the  width  of  the  plate,  the  load  bearing  area  will  be  set 
equal  to  the  plate  area.   Also,  at  this  point,  the  loading  pressure  will 
be  shifted  from  a  constant  ice  crushing  pressure  to  a  pressure  equal  to 
the  total  load  divided  by  the  plate  area. 

Initial  Assumptions  for  Rectangular  Plates 

The  plate  material  will  be  isotropic  and  homogeneous.   All 
four  sides,  of  the  plate  V7lll  have  a  clamped  edge  condition  for  support 
as  shown  in  Figure  2  and  only  in-plane  stresses  are  to  be  considered. 
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INTERACTION   AT 
PLASTIC  COLLAPSE 


Figure    1.      Idealized   Ice-Load 
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Figure  2.   Rectangular  Plate  Dimensions  and  Load  Position 
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Only  deflections  after  the  plastic  collapse  load  is  reached  will  be 
considered,   A  rigid,  perfectly  plastic  material  is  assumed  for  the 
collapse  mechanism  and  a  biaxial  strain  condition  will  be  assumed  at 
the  mid-point  of  the  plate. 

Model  of  Ice-Loads 

The  idealization  of  an  ice-load  will  be  a  four  sided  pyramid 
of  ice.   The  actual  applied  load  will  be  in  the  form  of  a  force  acting 
on  the  base  of  the  pyramid.   The  force  will  be  transmitted  to  the  rec- 
tangular plate  in  the  form  of  a  uniformly  distributed  load  acting  over 
an  area  centered  at  the  mid-point  of  the  plate.   The  value  of  the  load- 
ing pressure  will  be  the  ice  crushing  pressure  except  during  initial 
loading  and  after  the  bearing  surface  is  equal  to  the  plate  surface 
area.   These  exceptions  will  be  discussed  later. 

As  the  driving  force  increases,  the  bearing  surface  of  the 
pyramid  will  increase  as  the  ice  crushes.  Thus,  if  the  ice  crushing 
pressure  is  considered  constant  throughout  the  loading  process,  the 

bearing  surface  area  will  increase  directly  with  the  driving  force. 

2 
Since  the  area,  C  ,  at  any  time  is  equal  to  the  driving  force, 

F,  divided  by  the  ice  crushing  pressure,  P  ,  the  length  of  the  side  of 

the  bearing  area  will  be: 

This  relation  will  hold  after  the  pyramid  first  starts  to  crush  until 
the  plate  is  completely  covered  by  the  bearing  surface.   The  loading 
prior  to  the  initiation  of  ice  crushing  will  not  be  considered,  since 
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the  plate  deflection  will  be  in  the  elastic  range. 

Once  the  bearing  surface  equals  the  area  of  the  plate,  the 

applied  load,  F,  will  remain  constant  unless  P   is  allowed  to  increase. 

Thus  it  will  be  assumed  that  ice  crushing  ceases  and  P  will  increase 

c 

with  any  increase  in  F  as: 

P   =  F  (2) 

c 

AxB 

Plastic  Collapse  Load 

The  collapse  mechanism  for  the  plastic  analysis  is  the  general 
roof  shaped  deflection  surface.   For  the  transverse  cross  section  a 
three  hinge  collapse  mechanism  is  used.   For  any  loading,  it  was  assumed 
that  the  center  plastic  hinge  forms  either  simultaneously  with  the  two 
edge  hinges  or  at  some  increased  load  beyond  the  formation  of  the  edge 
hinges.  

Since  the  analysis  is  for  large  deflections  of  plates  the 
deflection  at  the  plastic  collapse  load  will  be  the  first  to  be  considered. 
Unfortunately,  there  is  no  easy  solution  for  the  deflection  at  this  load. 
The  problem  is  further  complicated  when  partial  loads  are  considered. 
Finite  displacement  relations  from  the  theory  of  plasticity  and  a 
membrane  solution  were  considered  to  be  the  easiest  methods  for  finding 
the  deflections,   A  similar  degree  of  accuracy  may  be  obtained  from 
either  method.   The  membrane  solution  has  been  shown  to  provide  good 
representation  of  the  plate  deformation  for  loads  well  beyond  the  collapse 
load.   Therefore,  it  was  decided  to  use  the  membrane  solution  throughout 
the  analysis. 
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As  the  pyramid  of  ice  is  forced  against  the  plate,  the  pressure 
acting  over  the  initial  bearing  surface  vjill  rise  until,  the  ice  crushing 
pressure  is  reached.   During  this  time,  the  plate  will  behave  elastically, 
From  the  time  P   is  first  reached  until  plastic  collapse  of  the  plate, 

the  area  of  the  bearing  surface  will  increase.   There  will  be  associated 

2 
with  the  collapse  load  F  ,  a  constant  P  ,  and  an  area  C   .   In  order  to 
^         o  c  o 

use  the  membrane  solution  to  find  the  deflection  at  collapse,  the  value 

of  C  must  be  knoiv'n.   Actually  it  is  easier  to  use  the  non-dimensional 
o 


form  C/A. 


The  moment  M   at  the  center  of  the  plate  was  found  using  an 

X 


elastic  solution  for  clamped  plates  under  partial  loads.   M  =  ^  P 

where  P  is  the  total  load  and  the  load  coefficient  (3  is  given  in 

2 

Figure  Al  for  values  of  B/A  and  C/A.   Then  by  setting  M  =  M  =(711  /4, 
*^  '^  X    o    o     ' 

the  yield  moment  for  a  rigid,  perfectly  plastic  material,  the  bearing 

surface  area  can  be  related  to  the  yield  stress  and  plate  thickness. 

2 
Thus,  when  P  =  P  C   and  both  sides  of  the  equation  are  divided  by  the 

width  of  the  plate,  A,  the  following  non-dimensional  equation  is 
2 


formed 


0-0  =  4  ^ 


Pc 


if 


(3) 


From  equation  (3)  the  bearing  surface  area  at  the  collapse 
load  may  be  determined  for  any  plate.   Various  solutions  to  the  equation 
are  given  in  Figures  A3-A5  for  given  values  of  B/A  and  A/H. 


2 

This  is  equation  (A7)  of  Appendix  A. 
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Plate  Deflection  Calculations 

Once  the  size  is  determined  for  the  load  bearing  surface  at 
the  collapse  load,  it  is  necessary  to  calculate  the  deflection  of  the 

plate  at  this  load.   As  discussed  previously,  a  membrane  solution  will 

3 

be  used  for  the  deflection  calculations.   The  membrane  formula, 

t    =  AqA^k  (4) 

rr3  s 

in  general  only  holds  for  a  uniformly  distributed  pressure  over  the 
entire  plate.   This  meant  that  some  geometric  relationships  were  re- 
quired in  order  to  use  the  membrane  formula  with  partially  loaded 

4 
plates.    It  was  assumed  that  the  plate  in  the  region  of  the  partial 

load  behaves  as  a  membrane  after  plastic  collapse  was  reached.   Then 

the  membrane  or  in-plane  forces  along  the  edge  of  the  loaded  area  are 

carried  through  to  the  edge  of  the  plate.   This  implies  that  the  slope 

of  the  deflection  profile  at  the  edge  of  the  loaded  area  is  the  same 

as  the  slope  at  the  edge  of  the  plate.   Thus  the  conical  shape  is 

represented  for  the  concentrated  loads  and  the  full  membrane  shape  is 

used  when  the  load  is  uniformly  distributed  over  the  plate.   It  is 

believed  that  this  provides  a  good  representation  of  the  true  deflection. 

The  deflections  for  the  concentrated  loads  may  be  less  than  they  should 

be,  but  the  accuracy  should  improve  as  the  load  bearing  area  increases. 

The  total  deflection  calculations  consist  of  three  basic 

steps.   First,  the  deflection  Z  for  the  loaded  portion  of  the  plate  is 

found  from  the  membrane  formula.   For  the  width  of  the  applied  load  C 


3 

See   Appendix  B. 

4 

See  Figure   3, 
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and  the  pressure  P  ,  Figure  3  shows  the  relation  of  Z  to  the  total 
deflection.   The  second  step  is  to  determine  the  slope  of  the  plate 
at  the  edge  of  the  loaded  area.   The  derivative  of  the  deflection 
profile  provides  the  slope  at  the  point  x  =  (A  t  C)/2.   Finally,  the 
total  deflection  at  the  mid-point  of  the  plate  is  found  by  combining 
the  membrane  solution  with  the  assumed  slope  at  the  plate  edge.   A 
simple  geometric  relation  will  determine  the  deflection  at  the  edge 
of  the  loaded  area,  and  the  membrane  solution  may  be  added  to  it  to 
give  the  total  deflection.   The  central  deflection  W  is  given  by: 

W  =  Z  -I-  2  (A-C)  tan  V   ,  (5) 

2 

where  the  membrane  formula 

Z  -  4  C^  k  Pp  (6) 

is  the  deflection  of  the  loaded  portion, 

and  tan  V     -  d_7^     =      P^  C  k^   at  x  =  A-C  ,  (7) 

dx         2S  2 

The  constants  k  and  k,  are  the  factors  that  depend  on  x.   The  values 

Q 

for  these  constants  are  included  in  Appendix  B  for  given  plate  aspect 

5 

ratios . 

Prior  to  discussion  of  the  procedure  for  load/deflection 
analysis,  the  membrane  solution  for  plate  deflection  must  be  understood, 
The  actual  formulation  of  the  deflection  relation  is  given  in  Appendix 
B.   There  are  two  dependent  variables  in  the  relation.   The  membrane 
force,  S,  V7hich  acts  along  the  edge  of  the  loaded  area  depends  on  the 
total  deflection  of  the  membrane,  W,  which  in  turn,  depends  on  S.   The 

See  Figure  Bl. 


usual  procedure  for  obtaining  the  value  of  S  is  to  assume  a  membrane 
deflection  and  shape.   Then  a  strain  relation  is  used  to  obtain  the 
principal  strain  in  the  plane  of  the  plate  for  the  assumed  deflection. 
The  usual  assumed  deflection  shapes  are  cylindrical  or  parabolic. 

In  order  to  model  the  ice  loading  condition,  the  equation  for 
the  shape  of  a  partially  loaded,  simply  supported  plate  was  used.   A 
biaxial  strain  relation  which  depends  on  the  B/A  ratio,  the  size  of  the 
bearing  surface  given  by  the  C/A  ratio,  as  well  as  the  maximum  deflection 
of  the  plate  was  established  for  the  mid-point  of  the  plate.    This  made 
it  possible  to  obtain  the  membrane  stress  from  the  stress-strain  curve 
for  the  plate  material.   Thus,  for  a  given  deflection,  the  membrane 
force  per  unit  length  is  equal  to  the  membrane  stress  times  the  plate 
thickness . 

If  the  membrane  form.ula  (equation  ''i)  is  used  alone,  the  applied 
load,  q,  could  be  found  that  would  produce  an  assumed  deflection.   Since, 
for  the  model  of  ice-loads,  the  loading  pressure  is  assumed  constant, 
an  iterative  procedure  must  be  used  to  obtain  the  proper  correlation 
between  Z  and  S  for  a  given  constant  load  P  . 

For  the  calculations  involving  the  deflection  of  a  partially 
loaded  plate,  the  total  deflection,  W,  is  used  instead  of  Z  to  determine 
the  value  of  the  membrane  force,  S.   This  appears  to  provide  a  more 
realistic  situation,  because  the  strain  relation  as  derived  in  Appendix 
D  depends  on  both  the  plate  shape  and  deflection. 


See  Appendix  D. 
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Load/Pef lection  Relations 

The  procedure  for  obtaiiiing  the  deflection,  W,  for  each  value 

of  the  total  load,  F,  will  be  iterative.   The  beginning  V7ill  be  based 

on  the  plastic  collapse  deflection.   After  the  value  of  (C/A)   is  found 

from  the  plastic  collapse  relation,  the  membrane  solution  will  be  used 

to  find  the  deflection  at  the  plastic  collapse  load,  F  .   The  first 

approximation  for  the  deflection  is  found  by  assuming  that  the  membrane 

stress  is  equal  to  the  yield  stress.   Then  for  the  value  found  for  C 

from  (C/A)   and  the  ice  crushing  pressure,  P  ,  the  first  values  of  Z, 
o  °    ^  '      c 

dZ/dx  and  W  (equations  6,  7  and  5)  can  be  calculated.   This  is  the 

starting  point  of  the  iteration  for  the  proper  value  of  W   for  the  given 

load  F  .   The  successive  solutions  of  W  will  follow  the  same  format  as 
o 

outlined  belo\i7  for  the  general  solution. 

For  any  given  total  load  F  ,  there  will  be  a  corresponding 

bearing  surface  area  represented  by  the  length  of  one  side  C  .   The 

problem  is  to  find  the  deflection  W.  that  corresponds  to  F   and  C  . 
'^  inn 

Let  F   =  F  +  A  F  where  the  subscript  n  designates  a  certain 

no  r  to 

value  of  static  load,  and  A F  is  the  incremental  value  of  the  total 
load.  The  subscript  o  refers  to  the  previous  value  of  applied  load. 
For  example,  F  will  initially  represent  the  collapse  load.   Then 

will  provide  the  necessary  load  bearing  area  to  determine  the  value 


of  W.. 

1 


The  previous  value  for  membrane  force  will  be  used  to  obtain 
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a  membrane  deflection  for  C  .   That  is, 

n 


Z.  -   ^  C^2   ^^  p^  ^5) 


rr3  s 

and  it  follows  that 


(i) 


Pc  Cn  kd  (10) 

2S 


and 

W.  =  Z.  -1-2  (A-Cn)  /  dZ  \  (11) 


?   =  Z   -1-  I    (A-Cn)  /^  dZ  \ 

2       \^  dx  y  i 


The  new  values  W.  and  C   are  used  to  enter  the  biaxial  strain  rela- 
1      n 


tions  and  obtain  a  new  value  for  the  membrane  force  S..   The  next  step 

1  ^ 

is  to  return  to  equations  (9),  (10),  and  (11)  to  calculate  Z.  -,  , 

(dZ/dx) .  , ,  and  W   ,,   The  subscript  i+1  is  used  to  distinguish  the  new 

values  of  deflection  from  the  membrane  force.   This  notation  indicates 

that  the  deflection  calculation  is  made  with  the  membrane  force  found 

from  the  previous  deflection,  W..   After  each  new  value  of  W.  is  found, 

1  1 

its  corresponding  membrane  force  S,  is  determined  from  the  strain 
relations  and  compared  with  the  previous  value  of  S,   Thus,  the  calcu- 
lation will  be  concluded  when  the  two  values  of  the  membrane  force 
compare  favorably. 

The  above  procedure  is  repeated  for  each  new  value  of  F^-^. 
When  the  load  is  reached  which  produces  a  load  bearing  area  equal  to 
the  area  of  the  plate,  that  is  when  C  -  A,  then  the  loading  pressure 

P  will  cease  to  be  constant  and  increase  with  the  total  load.   For 
c 

all  subsequent  increases  in  load  C/A  --  1.0  and  P^,  =  Fj^/(A  x  B) .   This 
assumption  of  continuous  loading  on  the  plate  is  just  an  extension  of 
the  model  of  the  ice  loading  beyond  the  point  where  the  pyramid  of  ice 
has  been  fully  crushed. 
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Energy  Absorbed  Relations 

Once  the  load /deflection  relation  is  found,  the  energy 
absorbed  by  the  plate  should  be  determined  at  each  value  of  F  in  order 
to  evaluate  the  possibility  of  plate  failure.   The  energy  of  deflection 
vill  be  the  greatest  at  the  edge  of  the  plate  except  for  cases  of  very 
small  load  bearing  surfaces.   For  the  plate  thicknesses  usually  con- 
sidered for  ships,  the  load  bearing  areas  will  be  large  enough  to  rule 
out  the  shear  effects  of  the  concentrated  loads. 

Since  the  plate  failure  can  be  expected  at  the  mid-point  of 
the  long  side  of  the  plate,  the  energy  absorbed  at  this  point  should  be 
greater  than  that  at  any  other  point.   There  will  be  two  sources  of 
strain  energy  at  the  edge  of  the  plate.   The  initial  source  will  be  from 
the  bending  moment.   Then  as  the  membrane  forces  increase,  the  portion 
of  the  total  strain  due  to  the  bending  moment  will  decrease.   The  com- 
bination of  the  strain  due  to  bending  and  stretching  will  give  the 
total  strain  at  the  expected  point  of  failure. 

Using  the  assumption  that  plane  cross-sections  remain  plane 
during  the  bending  and  stretching,  the  strain  distribution  across  the 
plate  thickness  will  be  as  shown  in  Figure  4.   The  rotation  of  the 
cross-section  is  due  to  the  plate  bending,  and  the  location  of  the 
neutral  axis  away  from  the  mid-plane  is  determined  by  the  size  of  the 
membrane  stress.   The  influence  of  the  membrane  stress  is  equivalent  to 
lowering  the  neutral  axis  from  the  geometric  center  to  a  new  point  z  -|^. 
This  shift  in  axis  of  rotation  of  the  cross-section  depends  on  the 
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Figure  4.   Strain  Distribution  for  Combined  Bending  and 
Membrane  Forces 
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ratio  of  the  membrane  stress  to  the  equivalem  yield  stress- 

H  _0 

2  a 


Zl  =  H  _^  •     .  (12) 


P 

where     O^   =  _      CT^  =   1.156   O  ^    for    Z/ =   0.5  (13) 

Vl-    V-    1^2 

This  relation  will  hold  true  for  values  of  0"^,  greater  than  0  since 
at  that  point  the  neutral  axis  will  lie  "outs  Me"  of  the  material  and 
all    strains    across    the    thickness   v;ill   be    in   tension. 

The    strain   of    the    outermost    fiber   al     the   mid-point    of    the 
long    side   v/ill   be    represented   by     6     .      From  Fi};ure   4 

62-65+21   tan    0  (]^4) 


Since  6   =  Zi  tan  0  =  H  (3^    tan 

=      '  2-ap 


6   -Htan©/    a„+l 

^2         i  -^ 


(15) 

Most  experiments  have  shca^/n  that  the  longitudinal  strain 
6   ,  at  the  mid-point  of  the  long  side  of  the  plate  is  very  small 
compared  to  the  transverse  strain  at  that  point.   Thus  £       was  taken 
as  the  principal  or  transverse  strain  at  the  point  of  expected  failure, 
Again  using  a  biaxial  strain  relation  as  given  in  Appendix  C   the 
uniaxial  strain  is, 

^^'-^,  (16) 

where  1/  is  equal  to  Poisson's  Ratio. 


7 

Young  in  reference  E5  recommended  the  use  o\'     (J      in  the  equation 

for  the  shift  of  the  neutral  axis  due  to  the  a.lJition  of  membrane 
stresses . 
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In  order  to  find  the  stress  at  the  expected  point  of  failure, 
the  uniaxial  strain  V7as  used  in  the  constituative  relation  for  the 
plate  material.   The  total  stress  will  be  designated  O      and  may  be 
found  for  the  strain  6  ^^  from  the  stress-strain  diagram  of  the  material 
in  simple  tension. 

The  energy  absorbed  per  unit  volume  at  the  mid-point  of  the 
long  side  will  then  be  represented  by: 


u  -    e  ,  0- ,    |. 


imibs  (17) 

in^   / 


This  value  may  be  found  at  any  point  in  the  loading  sequence  if  it  is 

assumed  that  6,  the  angle  of  rotation  of  the  cross-section,  is  equal 

to  y*  ,  the  angle  of  the  slope  of  the  deflection  profile  at  the  edge 

of  the  plate  as  shown  in  Figure  3.   This  is  a  reasonable,  as  well  as 

a  convenient  assumption,  since  the  slope  at  the  edge  is  determined 

during  the  deflection  calculations.   Thus,  for  each  value  of  F  and  its 

corresponding  W, 

tan  e  -  dZ  (18) 

.    dx 

Then   6  .  may  be  rewritten  as 

e.  =      H        /  dz  W  _o;s  +  i\  (19) 

2(1-1^^)  \  dx  /  \   0-p     / 

Also,  the  stress-strain  relations  for  the  plate  material  that  are 
required  for  the  deflection  relations  may  be  used  to  find  (J  . •   There- 
fore, for  any  loading  and  corresponding  deflection,  the  energy  absorbed 
by  the  plate  at  the  mid-point  of  the  long  side  may  be  determined.   A 
comparison  of  this  energy  to  the  restraint  energy  of  the  material  will 
lead  to  a  failure  criterion. 
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Failure  Criterion 

When  a  plate  is  drawn  in  biaxial  tension,  the  strain  rate 
becomes  unstable  when  the  principal  stress  reaches  the  ultimate  strength 

of  the  material.   The  plate  begins  to  decrease  in  thickness  rapidly  like 

(Pi) 

a  tension  specimen  begins  to  neck  at  the  ultimate  stress.       For  this 

reason  it  was  decided  to  use  the  ultimate  stress  of  the  plate  material 
to  establish  the  failure  criterion.   The  ultimate  stress  will  be  defined 
as  the  stress  at  v;hich  the  slope  of  an  engineering  stress-strain  diagram 
is  equal  zero.   The  area  under  the  stress-strain  curve  for  all  strains 
less  than  the  strain  at  the  ultimate  stress  will  represent  the  stability 
of  the  plate  material. 

A  restraint  energy  relation  will  be  defined  as 

where  O       is  the  ultimate  stress  and   G   is  the  ultimate  strain, 
u  u 

Both  values  are  defined  at  the  point  of  zero  slope  of  the  stress-strain 
curve . 

The  term  restraint  energy  is  used  since  R  represents  the 
amount  of  energy  that  might  be  absorbed  by  the  niaterial  before  insta- 
bility begins.   A  more  accurate  measure  of  the  ability  of  the  material 
to  absorb  energy  would  be  the  area  under  the  stress-strain  curve.   Since 
this  area  is  usually  difficult  to  obtain  readily,  the  R  relation  is  used, 
It  provides  a  reasonable  representation  of  the  area  under  most  stress- 
strain  curves.   Furthermore,  it  is  more  desirable  to  obtain  a  relation 
that  will  give  a  consistent  measure  of  the  restraint  energy  rather  than 
an  accurate  one. 
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The  two  parameters  that  best  correlate  the  minimum  amount 
of  failure  data  that  is  available  are  the  deflection  to  span  ratio, 
W/Aj  and  the  absorbed  energy  to  restraint  energy  relation  U/R.   These 
non-dimensional  parameters  represent  virtually  all  of  the  factors 
involved  in  the  plate  problem.   Several  other  parameters  were  considered, 
but  there  did  not  appear  to  be  any  relation  between  them  and  the  failure 
data.   The  W/A  ratio  helps  to  bring  together  all  the  geometric  proper- 
ties of  the  p].ate.   The  material  properties  are  well  represented  by  the 
U/R  ratio. 

Figure  5  is  a  plot  of  the  two  non-dimensional  parameters  for 
most  of  the  rectangular  plate  failure  data  that  is  available.   The  curves 
are  drawn  through  data  points  obtained  by  the  calculation  procedure 
described  in  the  previous  sections.   The  failure  points  for  each  plate 
are  designated  in  two  positions.   The  calculations  for  the  deflections 
of  the  plate  did  not  correspond  directly  with  the  experimental  data  as 
given  in  references  E2  and  E3 .   In  general,  the  calculated  deflections 
were  slightly  lower  than  the  experimental  deflections  for  a  given  total 
load.   The  possible  reasons  for  this  discrepancy  V7ill  be  discussed  later. 
It  was  felt  that  both  failure  deflection  and  failure  load  should  be 
indicated  in  Figure  5.   Thus,  the  first  set  of  failure  points  were  found 
by  using  the  experimental  deflection  at  the  mid-point  of  the  plate  to 
determine  the  W/A  ratio  at  failure.   The  second  set  of  points  were  ob- 
tained by  using  the  failure  load  from  the  experimental  data.   For  each 
rupture  pressure,  the  total  loading  force  F  was  calculated  and  from 
Figure  8,  the  corresponding  U  was  found.   The  U/R  values  for  the  failure 
loads  are  shown  on  Figure  5  as  the  lower  set  of  failure  points  on  the 
curves.  i 
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The  lower  set  of  failure  points  was  chosen  to  represent  the 
failure  data.   Although  the  lower  set  provides  the  more  conservative 
choice  for  a  failure  criterion,  the  major  reason  for  using  this  set  is 
that  the  failure  data  corresponds  better  to  the  deflection  and  absorbed 
energy  calculations  described  in  the  previous  sections. 

A  line  can  be  drawn  on  Figure  5  connecting  most  of  the  fail- 
ure points.   With  the  exception  of  the  plate  designated  IB  this  line 
represents  a  minimum  relation  of  VJ/A  to  U/R  for  plate  failure.   If  it 
is  moved  slightly  to  the  left,  the  failure  point  for  plate  IB  will  also 
lie  above  the  line.   This  new  position  of  the  limit  line  will  be  used 
as  the  failure  criterion. 

If  S  =  W/A  and  0  =  U/R,  then  failure  may  be  expected  beyond 
the  line  defined  by 

6  =  .278  -  0.182  0  (21) 

Greenspon  in  reference  T4  used  essentially  the  same  failure 
data  and  stated  that  failure  should  be  expected  for  loads  beyond  the 
pressure  that  produces  W/A  =»  0.10.   This  limit  as  shown  in  Figure  5  is 
quite  conservative  for  0  <  0.9,  but  for  plates  with  6  <  0. 10  and  0  >  0,9 
the  limit  may  be  too  large. 

The  energy  relation  was  chosen  for  this  analysis  in  order  to 
provide  a  parameter  that  will  accurately  represent  all  types  of  loading 
conditions.   Unfortunately,  the  available  failure  data  is  for  loads 
uniformly  distributed  over  the  entire  plate.   In  order  to  properly 
relate  the  failure  for  partial  loads  and  uniform  loads,  equivalent 
absorbed  energy  relations  are  used.   Thus  the  failure  criterion  should 
hold  for  all  but  concentrated  loads  with  small  bearing  surfaces. 


35 


CO 

m 

D 

z 

< 

OJ 

^ 

h 

CO' 
U. 

X 

Ll) 

< 

d: 

2 

> 

z 

O 

o 

o 

D 

cx 

G 

111 

Ll 

UJ 

(0 

CO 

< 

(f) 

< 

CO 

O 

LiJ 

^ 

Lii 

nr 

a. 

3 

UJ 

<: 

<r 

-I 

1 

^ 

Q. 

X 

o 

CO 

M 

-H 

to 

W 

(0 

a) 
u 

c 

0) 

<w 
<u 

*t-i 
o 

m 
<0 
■u 
(0 


c<3 

a 
I 

0) 
CI. 

0) 

■u 

^^ 
o 


CO 

(0 
0) 

■u 

(U 
u 

t-i 
•I^ 
cd 

0) 
4J 

(d 
ru 


IT) 

•H 
Em 


36 


Summary  of  Method  of  Analysis 

The  method  of  calculation  of  the  deflection  and  the  absorbed 
energy  for  any  plate  has  been  incorporated  into  a  short  computer  program. 
.The  description  of  the  program  itself  is  in  Appendix  E.   The  only  dif- 
ference between  the  procedure  summarized  belovv'  and  the  computer  program 
is  the  method  of  determining  a  stress  for  a  given  strain,   A  generalized 
equation  for  the  plastic  portion  of  the  stress-strain  curve  is  used  in 
the  program.   With  the  proper  choice  of  variables  the  generalized  stress- 
strain  equation  can  represent  the  actual  stress-strain  curve  of  a  given 
material  v/ith  reasonable  accuracy.   There  are  several  variations  to  the 
general  equation  that  may  be  used.   The  relation  described  in  Appendix 
E  appeared  to  be  one  of  the  easiest  to  use. 

The  computer  provided  the  speed  and  accuracy  to  carry  out 
the  calculations  for  many  plates.   The  hand  calculations  are  not  diffi- 
cult, but  they  are  quite  tedious.   With  the  use  of  the  tables  and  curves 
presented  in  the  Appendixes,  the  time  for  calculation  can  be  shortened 
considerably.   The  step-by-step  calculations  for  a  plate  of  knov^^n  geo- 
metric and  material  properties  and  a  given  ice-crushing  pressure  are: 

1,  For  the  values  of  B/A,  A/H,  and  C  /P  enter  the  curves 

o  c 

given  in  Figures  A3,  A4  or  A5  to  obtain  the  value  of  C/A  at  the  collapse 

2 
load.   Then  calculate  F   ^  C   P  . 

o       c 

2.  Calculate   Z   -   4   C^k   ?r.      ,    at'id      dZ    =      P^    C   kd    , 

n-3   S  dx  .2S 

and  W  =  Z  -f  JL  (A-C)   dZ  where  k  and  k^  may  be 
2        dx 

obtained  from  Figure  Bl  for  B/A.   Also  the  membrane  force  for  the  first 
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approximation  will  be  assumed  to  equal  the  yield  stress  tim.es  the 
thickness. 

3.  The  quantities  W/A,  B/A,  and  C/A  are  required  to  find 

the  value  of  the  average  transverse  strain  at  the  mid-point  of  the  plate. 
Tables  of  values  of  the  average  transverse  strain  and  representative 
curves  are  given  in  Appendix  D. 

4.  The  new  value  of  membrane  stress  will  be  found  from 
the  stress-strain  curve  for  the  plate  material  for  the  average  trans- 
verse strain  found  in  step  number  3.   Then  S,-  -    Cf      Yi. 

5.  Steps  2  J  3,  and  4  are  now  repeated  as  described  on  page 
28  until  the  new  value  of  S .  is  within  an  accepted  difference  from  the 
previous  value.   (An  absolute  difference  of  25  psi  is  used  as  a  tolerance 
in  the  computer  program) . 

6.  For  the  final  values  of  O's  a^id  dZ/dx  from  step  5,  cal- 
culate  e.  -     H     f    dZ    \    (      CTs  ^  -V 

7.  Again  enter  the  stress-strain  curve  for  the  plate  material 
to  find   0"^  which  corresponds  to   6^. 

8.  Calculate  the  energy  absorbed  at  the  point  of  expected 
failure.   U  =  6  j-  (T  t 

9.  Use  the  value  for   6  „  as  defined  on  page  33,  and  calcu- 
late R  -  6  ^   C^. 

10.   Using   S  =  W/A  and  0  =  U/R,  solve  the  failure  criterion 
equation  FC  -  S  +  .182  0  -  .278.   If  FC  >  0,  then  failure  can  be 
expected  for  any  additional  loading. 
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11.  If  FC  *==  0,  increase  the  total  load  F  by  the  load  incre- 
ment A  F  which  may  be  any  predetermined  amount.   Obtain  a  new  value 
for  C  -  Vf/P^, 

12.  Repeat  steps  2  through  11  to  find  the  deflection  and 
absorbed  energy  values  for  the  new  loading  condition. 

It  is  not  necessary  to  evaluate  the  possibility  of  failure 
at  each  load  increment.   After  one  or  two  times  through  the  procedure 
the  approximate  point  of  failure  may  be  apparent  and  the  solution 
method  can  be  shortened . 
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RESULTS  OF  APPLICATIONS  OF  THE  MEIBOD  OF  ANALYSIS 
Plastic  Collapse  Load  and  Deflection 

The  method  for  determining  the  total  load  for  plastic  collapse 
depends  on  a  given  ice  crushing  pressure.   The  area  over  which  this 
pressure  must  act  to  cause  collapse  is  determined  as  outlined  on  page  22, 
If  this  method  is  reversed  so  that  for  a  given  C/A  and  A/H  a  certain 

0"  /P^  is  found,  then  the  required  pressure  for  collapse  will  be  known. 
When  this  pressure  is  compared  to  elastic  experimental  data  for  plates 
of  the  same  size  and  material  properties,  the  collapse  pressure  should 
be  considerably  greater  than  the  elastic  loading  pressure.   It  is  very 
difficult  to  compare  deflections,  although  the  deflection  should  be 
increased  by  an  amount  similar  to  the  loading. 

The  plastic  collapse  load  and  deflection  were  calculated  for 
plates  of  the  same  size  and  material  properties  as  those  used  by  Sturm. 
and  Moore   ''   ,   The  size  of  the  bearing  surface  was  fixed  at  2"  x  2" 
and  12"  x  12".   Then  the  plastic  collapse  load  was  found  from  Figures 
A3-A5.   In  this  case  P  ,  which  has  been  called  the  ice  crushing  pres- 
sure, is  the  collapse  load.   Finally,  using  these  values  for  C  and  P^ 
the  collapse  deflection,  Wq,  was  determined  by  steps  2-5  of  the  previ- 
ous section.   The  results  of  the  calculations  for  three  of  the  plates 
are  shov7n  in  Table  1  along  with  experimental  data  from  reference  E4. 

A  second  comparison  of  the  plastic  collapse  loads  and  deflec- 
tions was  made  with  the  results  from  the  large  deflection  theory  of 
plates.   This  provided  a  means  of  judging  the  accuracy  of  the  membrane 
solution  method  for  calculating  the  plastic  collapse  deflections  for 


40 


uniformly  loaded  plates. 

Reference  E8  compares  many  of  the  solutions  to  the  large 

deflections  of  plates  to  experimental  results.   Most  of  the  relations 

are  given  in  the  form  of  curves  for  the  non-dimensional  deflection, 

"stress,  and  pressure  ratios.   The  deflections  of  three  plates  were 

calculated  from  the  corresponding  pressure  ratio.   The  plates  were  of 

the  same  size  and  material  as  those  given  in  Table  1.   The  collapse 

pressure  obtained  from  Figures  A3-A5  for  C/A  =1.0  was  used  in  the 

pressure  ratio, 

P.R.  =  Po^^  (22) 

T6EH^ 

For  example,  Figure  A3  shows  that   O'^/P^  =  215  for  an  aluminum 

(E  =  lOxlO^psi)  plate  I"x48"x48"  with  0      =   38800  psi  and  C/A  -  1.0. 

o 

Thus  P   =  Pq  =  180.5  psi.   Using  the  value  in  the  pressure  ratio,  a 
deflection  ratio  W/H  =  0,9  was  obtained  from  Figure  24  of  reference  E8. 
The  same  procedure  was  folloived  for  the  other  two  plates.   The  results 
from  the  large  deflection  theory  and  the  membrane  solution  for  collapse 
deflections  and  loads  are  shol^7n  in  Table  2. 

The  Membrane  Solution  for  Plate  Deflections 

It  was  relatively  easy  to  determine  the  accuracy  of  the  mem- 


brane solutions  for  the  plate  deflection.   The  plate  dimensions  and 

mate 
(El) 


(E2  E3) 
material  properties  of  the  test  plates  used  by  Day     '     and  Loeser 


were  put  into  the  computer  program  to  determine  the  load/deflection 
relationships.   The  results  of  the  calculations  as  well  as  the  input 
data  are  shown  in  Tables  El  to  Ell  of  Appendix  E.   It  should  be  noted 
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that  the  input  variable  for  ice  crushing  pressure  for  these  calculations 
is  only  the  loading  pressure  that  is  used  for  the  first . calculation. 
The  pressure  used  for  any  additional  loading  is  equal  to  the  applied 
force  divided  by  the  plate  area.   It  V7as  necessary  to  follow  this  pro- 
cedure because  all  of  the  test  data  is  for  uniformly  distributed  loads. 
Actually  the  solution  method  is  the  same  as  if  an  ice-load  is  used  and 
C/A  =  1.0.   In  this  way  the  computer  program  may  be  used  for  any  type 
of  lateral  load  except  highly  concentrated  loads. 

Also,  the  first  calculation  for  deflection  given  in  Tables 
El  to  Ell  does  not  represent  the  collapse  deflection.   The  starting 
point  used  for  these  calculations  was  chosen  for  convenience  in  com- 
paring the  results  to  the  respective  experimental  data.   For  example, 
the  rupture  pressure  given  in  Day's  experiments  was  divided  by  five  and 
the  result  was  used  as  the  input  for  ice  crushing  pressure  in  the  com- 
puter program.   Then  the  increment  of  applied  force   ^  F  was  set  equal 
to  the  ice  crushing  pressure  times  the  plate  area. 

The  results  of  the  calculations  are  shown  in  Figvires  6  and  7. 
The  experimental  results  of  Loeser  and  Day  are  also  plotted  in  each 
figure. 

Energy  Absorbed  Relation 

It  was  of  interest  to  plot  curves  of  total  load  F  to  the 
energy  absorbed  U  for  three  sets  of  data.   The  F  versus  U  relations 
for  the  calculations  on  Day's  test  plates  as  given  in  Tables  El  to  E7 
are  shown  in  Figure  8.   Similar  curves  for  calculations  on  Loeser's 
test  plates  (Tables  E8  to  Ell)  are  plotted  in  Figure  9.   The  results 
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of  calculations  made  on  arbitrary  plates  under  idealized  ice-loads 
are  given  in  Tables  E12  to  E19,  and  the  F  versus  U  curves  are  shown 
in  Figure  10. 

The  calculations  for  Loeser's  and  Day's  test  plates  have 
'been  previously  described.   The  results  shown  in  Tables  E12  to  E19 
are  for  four  plates  with  the  same  B/A  ratio,  but  of  varying  thickness. 
Also,  two  sets  of  material  properties  were  used  for  each  plate.   The 
calculations  follow  the  format  as  described  in  the  method  of  analysis 
and  Appendix  E. 

A  particular  case  is  encountered  in  the  calculations  for 
the  arbitrary  plates  that  was  bypassed  for  the  calculations  of  Loeser's 
and  Day's  test  plates.   The  assumption  that  the  ice  crushing  pressure 
changes  when  C/A  -  1.0  has  been  discussed  previously.   The  analysis  of 
the  arbitrary  plates  required  this  change.   This  can  be  seen  by  tracing 
do\;7n  the  C/A  column  in  Tables  E12  to  E19.   For  all  but  the  thicker 
plates  C/A  starts  at  the  value  for  (C/A)q  and  increases  until  C  =  A. 
For  square  plates  the  ice  crushing  pressure  will  change  to  the  total 
load  divided  by  the  plate  area.   But  for  plates  where  B/A  >  1.0,  this 
change  can  cause  a  decrease  in  the  crushing  pressure  rather  than  the 
desired  increase.   As  a  consequence  the  deflection  will  be  smaller  than 
the  previous  deflection.   In  order  to  help  eliminate  this  unnatural 
discontinuity,  it  was  assumed  that  the  loading  surface  continued  to 
grow  after  C/A  =  1.0.   That  is,  the  loading  pressure  was  kept  constant 
until  not  just  the  V7idth  of  the  plate  but  the  entire  plate  was  covered. 
Then  for  any  further  increases  in  the  load,  the  pressure  increased  as 
the  force  divided  by  the  area  of  the  plate.   This  scheme  would  be  the 
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natural  procedure  to  follov7  if  the  calculations  were  being  done  by 
hand  but  it  required  a  few  extra  statements  in  the  computer  program. 

Application  of  Failure  Criterion 

The  relation  between  $  ^  W/A  and  0  =  U/R  has  been  plotted 
in  Figure  5  for  the  results  shown  in  Tables  El  to  E7.   The  failure 
criterion  was  established  from  these  curves  and  the  failure  points 
given  by  Day     '  "  .      The  establishment  of  the  6  versus  0  curves  is 
discussed  on  pages  33-35. 

For  the  W  and  U  values  given  in  Tables  E8  to  Ell,  the  corres- 
ponding o's  and  0's  were  calculated.   The  relation  between  the  two 

failure  parameters  is  shoiv'n  in  Figure  11.   These  curves  represent  the 

^  1     1   ^  T       (El) 
experimental  work  of  Loeser 

Similarly,  Figure  12  is  a  plot  of  8  and  C5  for  the  calcula- 
tions made  on  the  arbitrary  plates  given  in  Tables  E12  to  E19.   The 
arbitrary  plates  by  their  sizes  and  properties  represent  typical  hull 
plating  on  ice-strengthened  ships.   Thus,  the  curves  of  Figure  12 
give  a  good  indication  of  the  possibilities  of  failure  of  these  plates 
under  the  idealized  ice-loads. 
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Figure  8.   Load  versus  Absorbed  Energy  for  Day's 
(E2,  E3)  Experimental  Test  Plates 


49 


600 


5O0 


(0    400 

CO 

_J 


? 

o 

X  300 
Li_ 


200 


lOO 


O- 


PLATE 

NOS.  FRO  171 

REF,  El 

ALL  PLATES  - 

A=13  IKI. 

PLATE 

H  (IN.) 

B/A 

Z 

0.1  56 

a.s 

3 

0.1  SS 

1.5 

5 

0.166 

LO 

4 

0.^50 

2.0 

3  4  5  6 

Uxio-^GW) 


6 


10 


Figure  9.   Load  versus  Absorbed  Energy  for  Loeser's 
(El)  Experimental  Test  Plate 
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Figure  10.   Load  versus  Absorbed  Energy  for  Arbitrary 
Plates  of  Tables  E12  to  E19 
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DISCUSSION  OF  RESULTS 
Load/Pef lection  Relations 

Figures  5  to  11  and  Tables  1  and  2  are  presented  so  that  a 
comparison  of  the  various  portions  of  the  method  of  analysis  may  be 
made  with  associated  experimental  data.   The  degree  of  accuracy  of  each 
step  of  the  analysis  is  shown  in  the  tables  and  figures. 

Although  the  data  presented  in  reference  E4  V7as  published  in 
1937,  it  proved  to  be  the  most  useful  experimental  data,  on  rectangular 
plates  under  partial  loads.   The  comparison  between  Sturm  and  Moore's 
experiments  and  the  calculated  values  for  collapse  deflection  and  load 
given  in  Table  1  is  very  crude.   The  experimental  data  is  for  the 
elastic  range  only,  and  the  calculated  values  assume  that  the  plastic 
collapse  load  has  been  reached.   Hotvever,  the  experimental  data  does 
provide  at  least  a  lower  bound  to  the  calculated  values. 

There  are  two  general  points  that  should  be  recognized  from 
Table  1.   The  relation  between  the  deflections  and  the  B/A  ratio  that 
is  indicated  by  the  experimental  data  holds  for  the  calculated  values. 
That  is  J  the  experimental  deflections  increase  as  B/A  increases  and 
the  same  holds  true  for  the  calculated  values.   A  similar  comparison 
between  deflections  and  plate  thicknesses  shows  that  the  experimental 
deflections  increase  as  the  plate  thicknesses  decrease.   However,  the 
calculated  deflections  do  not  follow  the  same  trend  as  the  experimental 
results.   It  seems  that  the  thinner  plate  should  have  the  greatest 
deflection  as  it  does  in  the  experimental  data.   Although  only  relative 
comparisons  have  been  made,  the  failure  of  the  membrane  solutions  for 


54 


the  collapse  deflections  to  compare  favorably  with  the  experimental 
results  should  be  considered  a  limitation  of  the  method  of  analysis. 

A  comparison  of  the  plastic  collapse  deflection  obtained  from 
the  membrane  solution  with  the  deflections  calculated  from  large  deflec- 
tion theory  may  be  made  from  Table  2.   In  general,  the  membrane  solution 
deflections  are  lower  than  those  obtained  with  the  same  pressure  from 
large  deflection  theory.   For  this  case  of  uniformly  distributed  loads, 
both  sets  of  calculations  shav  that  the  deflection  increases  as  the  B/A 
decreases.   Also,  the  deflection  for  the  large  deflection  theory  remains 
constant  for  the  three  plates  considered.   Again  the  relation  between 
deflections  and  plate  thicknesses  does  not  compare  favorably.   The 
membrane  solution  maintains  approximately  the  same  W/H  ratio  at  collapse 
for  all  thicknesses.   If  the  results  from  the  large  deflection  theory 
are  used  to  define  the  plastic  collapse  condition,  then  the  limitation 
in  the  method  of  analysis  as  discussed  for  Table  1  is  again  encountered. 

The  accuracy  of  the  membrane  solution  may  be  easily  seen  in 
Figures  6  and  7.   The  calculated  deflections  are  either  completely  above 
the  experimental  deflections  or  completely  below.   In  all  cases,  if  the 
plate  has  strength  properties  (see  Tables  El  to  11)  greater  than  mild 
steel  (  0'   =  38000  psi)  the  deflections  lie  below  the  experimental 
data.   For  the  higher  strength  plates,  the  deflections  are  all  below 
the  experimental  results.   The  one  exception  is  Day's  test  plate  2A  in 
Figure  7.   It  has  a  yield  strength  slightly  above  the  mild  steel  plates 
and  the  curves  for  the  calculated  and  experimental  deflections  coincide. 

Since  the  shapes  of  the  calculated  and  experimental  deflection 
curves  agree  very  well,  the  only  major  discrepancy  is  the  difference  in 
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the  deflections  as  noted  above.   The  reason  for  the  differences  is  the 
solution  method  for  obtaining  the  membrane  force  for  the  membrane  for- 
mula.  This  same  reason  may  be  cited  as  the  cause  for  giving  lower 
values  of  deflection  for  high  strength  materials  and  higher  deflections 
for  lower  strength  materials.   The  generalized  stress-strain  curve  used 
in  the  computer  program  usually  gives  values  of  stress  that  are  lower 
than  the  actual  value  for  the  same  strain.   This  will  in  turn  produce 
larger  membrane  deflections  since  the  niembrane  force  is  S  =  ^s^' 
Furthermore,  for  two  plates  of  approximately  the  same  thickness,  the 
membrane  stress  will  be  greater  for  a  plate  made  of  high  strength  metal 
than  one  with  a  lo\^7er  yield  strength.   Consequently,  the  membrane  de- 
flections for  the  high  strength  material  will  be  less  than  those  of 
materials  with  a  lower  yield  stress. 

For  the  general  discussion  of  the  application  of  the  method 
of  analysis,  Figures  13  and  14  illustrate  what  the  load/deflection 
relationship  is  for  hull  plating  under  the  idealized  ice-loads.   The 
curves  and  the  respective  results  of  the  computer  calculations  given 
in  Tables  E12  to  E19  show  that  for  small  load  bearing  surfaces  the 
collapse  deflections  appear  to  be  excessive.   The  collapse  deflections 
for  the  large  bearing  surfaces  are  usually  smaller  than  what  might  be 
expected.   For  the  uniformly  loaded  plates,  the  collapse  deflection  is 
not  necessary  for  the  load  deflection  analysis.   Thus,  for  plates  that 
require  a  collapse  load  greater  than  that  given  by  step  1  of  the  method 
of  analysis,  the  first  value  for  deflection  given  in  Tables  El5,  El8, 
and  E19  does  not  represent  the  collapse  deflection. 

For  the  partial  loading  conditions  there  is  no  direct  relation 
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between  the  deflection  at  one  value  of  the  applied  force  and  the  next. 
Thus,  once  the  uniform  load  is  reached,  the  deflection  should  be  within 
the  same  accuracy  as  the  plates  sha^;n  in  Figures  6  and  7.   This  fact 
alone  should  help  to  substantiate  the  loading  process.   In  all  cases, 
the  deflect  ion/ load  relation,  in  general,  follo\\'s  a  continuous  curve. 
There  are  exceptions  to  this,  but  they  are  caused  by  the  change  over 
from  the  constant  ice  crushing  pressure  to  the  uniformly  distributed 
load.   If,  for  example,  the  deflection  for  the  first  partial  load  is 
greater  than  the  deflections  for  subsequent  loadings,  the  deflection 
will  be  decreasing  as  the  total  load  increases.   As  can  be  seen  in 
Figures  13  and  14,  the  deflection  is  never  less  than  the  previous  value. 
Therefore,  if  the  accuracy  of  the  calculated  deflections  for  the  uni- 
formly distributed  loads  can  be  accepted,  then  the  calculated  deflec- 
tions for  the  partial  loadings  should  be  of  reasonable  values  for  use 
in  the  failure  analysis  of  the  plates. 

Absorbed  Energy  Relation  and  Failure  Criterion 

There  is  one  basic  purpose  for  calculating  the  energy  absorbed 
by  the  plate  during  the  loading  operation.   The  calculated  energy  pro- 
vides a  means  of  determining  the  possibility  of  plate  failure  on  any 
increase  in  the  load.   There  are  several  methods  of  determining  the 
amount  of  energy  absorbed  by  a  plate.   The  method  presented  provides  a 
measure  at  the  point  of  expected  failure. 

It  has  been  found  in  experimental  tests  of  dynamically  loaded 
clamped  plates  (see  references  E19  and  E20)  that  a  linear  relation  exists 
between  the  measure  of  appi ied  force  and  the  measure  of  the  energy 
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absorbed  by  the  plate.   The  relation  between  F  and  U  in  Figures  8  and 
9  is  clearly  linear  for  all  plates.   The  results  shown  in  Figure  10  for 
the  arbitrary  plates  under  the  idealized  ice-loads  are  not  linear.   The 
nonlinearity  is  due  to  the  changing  size  of  the  bearing  surface  whereas 
the  curves  of  Figures  8  and  9  are  for  a  uniformly  distributed  load 
throughout  the  loading  process. 

It  is  interesting  to  note  that  the  only  difference  between 
the  two  curves  in  Figure  10  is  the  material  properties  of  the  plates. 
The  F  versus  U  curve  for  each  of  the  four  plate  thicknesses  considered 
coincides  with  the  other  plates  of  the  same  material  properties.   This 
same  general  relation  appears  to  hold  for  the  plates  shown  in  Figures 
8  and  9,  although  B/A  ratios  and  plate  sizes  differ  in  these  figures. 
Thus,  based  on  the  results  of  Figures  8  and  9,  and  the  consistency  of 
the  results  shown  in  Figure  10,  the  absorbed  energy  relation,  U,  should 
provide  a  good  measure  of  the  possible  failure  of  a  plate. 

The  calculated  values  for  8  and  0  from  the  data  given  in 
Tables  E8  to  Ell  are  plotted  in  Figure  11.   The  values  of  U  representing 

the  maximum  load  that  was  applied  to  each  of  the  experimental  plates  by 

/■pi  \ 
Loeser      are  also  indicated.   With  the  exception  of  plate  number 

five,  all  of  the  maximum  load  points  lie  on  or  above  the  failure  cri- 
terion line.   Plate  number  five  failed  due  to  a  weld  defect,  but  the 
other  plates  did  not  fail  in  Loeser 's  experiments.   Thus,  there  appears 
to  be  at  least  a  rough  correlation  of  the  failure  criterion  for  two 
different  experimental  results. 

Figure  12  shotjs  the  final  application  of  the  method  of  analysis, 
If  the  line  of  the  criterion  is  accepted  for  the  smaller  values  of   5 
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then  the  thicker  plates  can  be  expected  to  fail  shortly  after  their 
collapse  deflection  is  reached.   This  is  actually  what  would  be 
expected  for  the  very  stiff  plates.   However,  it  seems  that  the  line 
of  the  failure  criterion  should  include  larger  values  of  0  for  the 
smaller  values  of  O  .   This  is  just  a  speculation,  and  the  presently 
defined  failure  criterion  appears  to  provide  a  reasonable  estimate  of 
the  failure  load  of  a  plate, 

CONCLUSIONS 

The  method  of  analysis  for  large  deflections  of  a  clamped 
rectangular  plate  under  idealized  ice-loads  provides  a  reasonable  means 
of  predicting  the  behavior  of  a  plate.   Furthermore,  the  solution 
method  provides  a  simultaneous  check  on  the  possibility  of  plate  fail- 
ure.  The  calculations  are  tedious,  but  not  difficult,  and  the  use  of 
the  computer  program  will  greatly  aid  any  extensive  failure  analysis. 

A  failure  criterion  was  sought  and  obtained.   Furthermore, 
the  criterion  is  not  limited  to  plates  under  ice-loads.   It  should  be 
equally  applicable  to  the  analysis  of  a  plate  under  any  lateral  load 
except  in  cases  where  shear  effects  are  large. 

RECOMMENDATIONS 

The  major  recomrriendation  is  quite  obvious,  that  is,  to  conduct 
tests  to  provide  more  data  for  the  failure  criterioii.   Other  areas  that 
should  be  investigated  are  other  possible  solution  methods  for  the 
load/deflection  relations  and  the  energy  absorbed  relation.   These  areas 
should  be  included  in  any  failure  tests  since  they  are  necessary  for 
the  establishment  of  the  failuro  criterion. 
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Appendix  A 
PLASTIC  COLLAPSE  LOAD 

The  elastic  nioment  at  the  center  of  a  partially  loaded  simply 
supported  rectangular  plate  is  given  by  Tiinoshenko  in  reference  Tl.   The 
moment  about  the  longitudinal  axis  at  the  mid-point  is, 

M^  =   (3  P  (Al) 

where  P  -  qC^  and  0    is  given  in  Tables  20-22  of  reference  Tl.   The 
values  of  3   depend  on  C/A  and  B/A,  and  those  values  given  in  Tables 
20-22  decrease  in  accuracy  as  C/A  decreases.   Values  for  (3   for  the 
smaller  values  of  C/A  can  be  found  from  the  equation  for  M  given  on 

X 


p.  160  of  reference  Tl.   That  equation  is: 


M   =   P 

X 


8n 


(2  log  4A +  X  -  1.571)  (1+V)  +  p.   (l-V) 


(A2) 


The  parameters  X   and  /X.  depend  on  B/A  and  the  values  for  them  are 
given  in  Table  27  on  p.  161  of  reference  Tl.   The  values  of  |3   calcu- 
lated from  equation  (A2)  for  C/A  less  than  0.20  are  given  in  Table  Al. 

In  order  to  use  equation  (Al)  for  a  plate  with  a  clamped  edge 
condition,  the  coefficient  ^  must  be  adjusted,  since  the  values  given 
in  Tables  20-22  of  reference  Tl  and  Table  Al  are  for  a  simply  supported 
plate.  The  correction  term  Pi  may  be  found  on  p.  206  of  reference  Tl 
for  various  B/A  ratios.  Thus,  the  maximum  moment  at  the  center  of  a 
clamped  rectangular  plate  under  a  partial  load  is, 

(M  )  max  =  (  p  +  §i)   P  =   (3'P  (A3) 

The  factor   ^'  is  plotted  against  C/A  for  three  values  of  B/A 
in  Figure  Al.   Tv/o  steps  were  necessary  in  finding  the  proper  value 
of  p   .   First,  P   was  plotted  against  C/A  for  the  values  of  ^   given 


in  Tables  20-22  of  reference  Tl  and  for  those  values  of  (3   given  in 
Table  Al.   Then  for  each  B/A,  a  single  curve  was  formed  which  combined 
the  two  solutions  of  0   .An  example  of  the  graphical  solution 
for   0   is  shown  in  Figure  A2 .   The  second  step  required  adding  the 
constant   3  x  to  the  new  values  of   ^  obtained  graphically  to  complete 
the  solution  of   ^ 

Once  the  elastic  solution  for  the  moment  at  the  mid-point  of 
a  clamped  rectangular  plate  was  established,  the  plastic  collapse  load 
was  found  by  simple  plastic  analysis.   The  plastic  moment  at  a  cross- 
section  through  the  thickness  of  a  rigid,  perfectly  plastic  plate  is 
Mo  =  Q'o  H^ 


(A4) 


Setting  (M  )  max  =  M^  the  plastic  collapse  load  becomes, 

Po  =  Oo   H^ 
4  P' 

2      2 
Since  Pq  =  qoC   -  P^^C   for  an  ice  crushing  pressure  equal  to  P^ 

C^  -  (To  H^ 


^^Pc 


or 


c   ~-  H  /g^ 

By  dividing  through  by  the  width  of  the  plate,  A,  equation  (A6)  may 
be  non-dimensionalized . 
C 


or 


A 
Co 


2A  Y^'Pc 


4  ^' 


A 
H 


(A5) 


(A6) 


(A7) 


The  desired  relation  between  the  load  bearing  surface  and  the  ratio  of 
material  and  ice  properties  for  given  plate  geometry  may  be  obtained  from 
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equation  (A7). 

The  parameter  Of^/P^  for  values  of  B/A,  C/A,  and  A/H  has  been 
tabulated  in  Tables  A2,  A3,  and  A4.   From  these  values  the  curves  of 
A/H  versus  Cq/P^  for  various  values  of  C/A  have  been  drawn  in  Figures 
A3~A5.   It  was  found  that  this  was  the  best  way  to  present  the  data  in 
order  to  make  use  of  them  in  the  plate  analysis.   For  a  rectangular 
plate  of  given  dimensions  and  yield  stress  being  loaded  by  ice  with  a 
given  crushing  pressure,  (C/A)^,  the  size  of  the  load  bearing  surface 
at  the  plastic  collapse  load  may  be  found  for  A/H  and  CT  /P^  from  the 
curves  of  Figures  A3-A5. 
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RATTH    GP    YlFLn    ST^^FSS    TH    ICF    CRUSf^MG    PKFSSIPF 

R/A^  loOOO 


A /I! 

10,0 

20  eO 

30,0 

40o0 

50^0 

lOOeO 

C/A 

OoO?0 

0,0  ■ 

Oo2 

0o4 

0.^ 

1.  ? 

4,9 

OeOAO 

0«? 

Oe6 

K4 

2o  6 

4,0 

16,1 

0,  060 

0.3 

1,3 

3.0 

5,3 

8,2 

3  2,  "^ 

0.080 

Oo  5 

2.1 

Ao^ 

806 

13,4 

5  3,  6 

0.  ICO 

0.  R 

?o  1 

7,1 

12,6 

19,  6 

7  8,6 

0.200 

2.5 

9.Q 

2?,^ 

39,7 

62o  0 

248,0 

OoAOO 

6, a 

27',  1 

6  ]  ,  1 

lOP,  5 

1  6  '3  ,  6 

67  8,4 

0,600 

10.  4 

^  1  •  .« 

^4,0 

167.0 

26]  .0 

lC4^o0 

0.800 

llo  5 

46,  ] 

103=  7 

184,3 

2  8  8,0 

1152o0 

1.000 

9,2 

3-^.0 

^3,2 

147,8 

2  31.0 

^24,0 

Table  A2 
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RATIO    TF    YIFLH     STRES"     T^     TCF     CRUSHING    f^RFSS'.J^.E 

P/A=  lo400 


A/H 

10.0 

20.0 

3  0.0 

40  eO 

50,0 

lOOoO 

C/A 

■o,o?o 

0.  1 

0.2 

0.5 

OoR 

]o3 

5o2 

0,0^0 

0.2 

0,7 

1.6 

2.^ 

4«3 

17,4 

0,060 

0.4 

1<.4 

3,2 

5.7 

8,0 

35„7 

0«0^0 

0„6 

2.4 

5.4 

Q.  f^ 

15.0 

60,0 

0.  100 

0.  q 

3,5 

P.O 

14.2 

22.  1 

8  8.  5 

0,200 

2o  R 

lU? 

25.2 

44,  ^■ 

70,0 

28O0O 

0.400 

7.8 

3  1  o  4 

7  0o6 

]  2  '^  .  4 

196.  0 

^84.0 

0.^00 

12,  A 

4^.5 

in,5 

1^8,  1 

300.  6 

12^8.4 

0,  noo 

14,8 

so,  1 

13  2,9 

2  3  6o  ^ 

^6^,3 

1477, 1 

1,000 

14.  0 

ss.  p 

12^06 

2  2  3  .  4 

^49.0 

1 396.0 

Table   A3 
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RAin    fJF     YT[-LD    STRCSS    TO     TCE     r^iJ^Hir'G    PRESSURE 

B/A=  .?oC00 


A/H 

10.^ 

20.^ 

30.0 

40,0 

5  0.0 

100,0 

C/A 

0,0?0 

Oo     1 

Oo2 

0,^ 

0,8 

lo? 

5,3 

0.04  0 

0.^ 

0,7 

1.6 

2.8 

4,4 

17,6 

0,060 

Oo4 

U5 

3.3 

5o  ^ 

9,  1 

3  6o6 

O.OPO 

0«6 

2o4 

5.5 

9,7 

15,2 

6  0.7 

0.100 

0.^ 

3,6 

8.  1 

14.3 

2  2,4 

P9o6 

Oo?CO 

2.Q 

11  o5 

^5.0 

^'6o  1 

7  2,0 

2  8  "  »  0 

0.400 

8.  2 

32,8 

7^,7 

1  31,  1 

2  0^c  8 

81^.2 

0.600 

13,4 

S3,  6 

12  0.5 

214.^ 

3  3  4o  R 

1 33^,2 

0,  ROO 

16„6 

66o6 

1-^^^.  S 

2  6  6,2 

^  16.0 

16  64,0 

1,000 

1  6  o  4 

65,6 

1  A7,  ^ 

?6^,4 

410,0 

164  0,0 

Table  A4 
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Appendix  B 

DERIVATION  OF  MEMBRANE  FORMULA  FOR  PLATE  DEFLECTION 

The  solution  for  the  deflection  of  a  uniformly  loaded  mem- 
brane is  given  in  several  references.   (See  for  example,  references 

Tl,  T4,  and  El).   The  method  presented  here  is  the  one  given  by 

(El) 
Loeser.       Since  the  method  has  been  developed  elsewhere,  only  the 

major  points  of  the  solution  are  given. 

The  differential  equation  of  a  uniformly  loaded  rectangular 

membrane    is: 

a^z  a^z         ^      _gL  (Bl) 

ex^     "^   a  y^        "    s 

vhere    q    is    the    lateral    loading    pressure   arid   S    is    the   membrane    force. 

For   the    plate   dimensions    given    in  Figure    2,    the   boundary   conditions    are: 

Z-OatX^O,    X=A,    y=+     B  

2 

The  assumed  solution  for  equation  (Bl)  is  the  series: 

oo 


Z  =    >         b^  sin   nn:x    Y^  (B2) 

£- /  A 

n«l,3,5 

The   b    's    are    constant    coefficients,    and   Y    's    are   dependent   upon  y    only, 
After  differentiating    tv/ice  with   respect    to  x   and   y   and    substituting 
f°^     9    Z      and     8^Z      ,    equation    (Bl)    becomes: 


oo 


n=l,3,5 


bn   si"  ilOjl-  f    _      "^^^      Y^  +  <     )  -    -  ^  (B3) 

a2  is 
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Expansion  of  the  right  hand  side  of  equation  (B3)  into  a 
Fourier  series  and  using  constant  value  for  q  between  0  <  x  <  A 
results  in  the  following  expression: 

oo  n-1 


3.   =   Aq     >         (-1)     sin  nrtx  (P4) 

S    rtS    ^ ^  „  _      n        A      , 

n  =  1,3  ,5 

Setting  equations  (B3)  and  (B4)  equal  to  each  other  gives 

n-1 

A^  S  nn-b„ 

n 

Solution  of  this  ordinary  differential  equation  for  the  conditions 
that  Y   is  symmetric  about  y  =  0  and  that  Y   =  0  when  y  =  +  B/2  yields 

n-1 

Y^    ^   4qA^    (-1)^ /    cosh   nTt  B        cosh   nfTy       ]         (B6) 

^   ^      3  ^3  ,       cosh  n/TB  2A~  A      ~      } 

S   n     a      b^  2A-      \  / 

Then   the   membrane   deflection  may  be    obtained   by    substituting   equation 

(B6)    into  equation    (B2). 
oo 

V"  n-1  ,    ^ 

Z   =      ^^-r-^    o    q      4qA^    (-1)  (  A         ]    sin  nTtx      (B7) 

""   ~    ^'-^'^        ^      3^3       "   I     ^"    cosh  n/^B     I  a 

S   n    IT  y  -y-  / 

If    the    constant    terms    are    removed   and    the    series    is    defined   as   k,    the 
membrane   deflection  equation   is: 

Z   =  4cjA^       k  (B8) 


sn^ 


and  c^ 


Z 


n-1 


cosh   n  rt  y 


k   =       ^^7-r    ^    .       (-1)^  sin  nUx       I    ,      A_     \  (B9) 

n-i,J,D  ^^  ^  ^         ^^g^  nrrB 
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The  values  for  k  have  been  plotted  in  Figure  Bl  for  three 
values  of  x/A.   The  use  of  i  (lose  curves  eases  the  burden  of  computation 
for  the  deflection.   For  ex-uiiple,  the  value  of  the  central  deflection, 
W,  of  a  uniformly  loaded  pLiie  will  require  the  value  of  k  correspond- 
ing to  the  B/A  ratio  of  the  plate  and  x/A  =  0.50.   This  same  value  of 
k  will  hold  for  all  values  oi  q  and  S. 

The  slope  of  the  m-'mbrane  surface  is  obtained  by  differen- 
tiating equation  (B8)  with  i  .-spect  to  x.   Thus,  the  slope  at  the  mid- 
point of  the  long  side  of  tin-  plate  is: 


oo 


cosh  n rr 


1 


A?L  ^   AaA_   /         (-1)      ,   A_  I  cos  nrcx    (BIO) 

9x    Sn2   ^Tl^3^5    ^2    \  "  cosh  nrfB 


2A 


Then  for  y  =  0  and  x  =  0  or  x  =  A 


4q.A 

"  srr'-' 

oo 

dZ 

ax 

r,     -     1 

-(Ell) 


Noting  that   ^T~Y' 3  ^     ly   ^    JJZ.   »  equation  (Bll)  becomes 
"     '  '    n^     8 


00 


or 


8x    2S      srt-'   /r-t  -^    s    n2  ^osh  nHB 


6Z  ^  5A  k^    at  K  =«  y  =  0  (Bl2) 

ax      2s 


where 


00 


•^d  =  1. 


rt2  ,  ^2   cosh   nrrJB 

'-^'^  2A 


The   coefficient    k^    is    also  ^tveu    in  Figure   Bl    for  values    of   B/A. 
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Appendix  i; 
BIAXIAL  STRAIN  I'KIATION 


Assume  that  a  given  element  ol  the  plate  has  only  normal  f 
acting  on  each  face. 


orces 


Sn-* 


For    the   case    of   uniaxial    loadinjr    let    S      =    S      =0        Then   the 

y    2     ' 


strain  components  along  each  axis  will  b.- 

^xl"^^l'   ^yl  ~"   ""^^xl'  ^'zi  =   ~"^^xl 


(CI) 


where  G^  represents  the  principal  strain  ^n  the  x  direction.   The 
strain  components  in  the  y  and  z  directinug  are  obtained  by  considering 


S   =  S   =  0  and  S  -   S      =0  respectively.   Thus, 
X    z  X    y 


£y2  "  ^2'  e^2 


y2   ^'z2 


V  £ 


y2 


and     6^3-  63,   6^3=   -^6^3,  (.^3-^   ~l/6^3 


(C2) 
(C3) 


If  all  the  forces  are  acting,  1  i,a  total  strain  in  the  x  direc- 
tion will  be, 


X 


X 


£1  -V(  6^2+  6,3) 


y2   -1^6,3 


(C4) 


For  plates  it  is  generally  assmn.jd  that  6   is  small  compared  to  the 

other  strains.   The  total  strain  in  the  x  direction  for  €   =0  will  be 

z 

^x  --^  ^1   -^^y2  (C5) 

Similarly,  the  total  strain  in  Lhe  y  direction  for  G   =0  will  be, 

2 

^y  ^  ^2   "■^'^xl  CC6) 

For  the  case  v^here  the  in-planr-  strain  is  uniaxial  in  the  x  direction 

e    =  0  and ,  hence ,  6   =1^e,  =1/6,.   From  equation  (C2)  G  0=^0 
y  .^xii  y22 

=  Ve^.      Substituting  into  equai.ion  (C5)  £   =  (    1-  V  ^)  G    .    (C7) 

This  relation  will  be  useful  in  obtaining  the  stress  at  a  given  point 
where  the  strain  in  the  x  direcLion  is  known,  and  the  strain  in  the  y 
direction  is  zero.   The  stress  .an  be  found  from  a  stress-strain  diagram 
for  simple  tension  by  using 

€1-   -^ 

ITv^  -  (C8) 

A  more  general  relatj,ui  is  required  when  there  is  strain  in 

both  the  X  and  y  directions.   'J  1  £   ^^o,  the  strain  in  the  x  direction 

y   • 

may  be  found  by  combining  equal  ions  (C5)  and  (C6) . 

-  ^i    -  ^  (  e^,    - 1/61) 

6^  =  (1-  V^)  G  1  -VGy  (C9) 

When  used  in  the  form 

€.  ^    e  X  +  V  €y 

1  -V^  (CIO) 

the    stress    in   the   x   direction  niiy  be    found    from  a   stress-strain   diagram 

for   simple    tension. 
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Appendix.  D 
PLATE  DEFLECTION  PROFILE  AND  AVERAGE  TRANSVERSE  STRAIN 

In  order  to  use  the  membrane  formula  for  the  calculation  of 
plate  deflections,  a  relation  betv7een  the  deflection  and  the  membrane 
force  must  be  established.   This  is  accomplished  by  assuming  a  plate 
deflection  profile.   Then,  for  a  given  deflection,  the  average  strain 
within  the  plate  can  be  calculated.   From  the  stress-strain  curve  for 
the  average  strain,  a  value  for  the  in-plane  or  membrane  stress  can  be 
found  which  is  converted  to  the  membrane  force  by  multiplying  by  the 
plate  thickness. 

For  this  analysis,  it  was  desired  to  use  an  assumed  deflec- 
tion profile  that  would  represent  the  various  loading  conditions. 

Thus,  the  solution  for  a  partially  loaded  simply  supported  rectangular 

(Tl) 
plate  as  derived  by  Timoshenko      vas  selected.   The  general  expres- 
sion from  page  138  of  reference  Tl  for  the  deflection  at  any  point 
along  the  x  axis  (see  Figure  2)  is: 
oo  m-_l 

cosh(a^-2a"^0 


w 


4£A. 


(-1)    sin  mnc    \    ^ 1 


^n^     m  -  1,3,5 


m' 


2A 


cosh  (X 


m 


m  '^   m^ 


—     ^   m      m       m 


m 


2cosh  OC 


m 


sin  m  rtx 


(Dl) 


mTTB     ,   -,    mrtC 
where  CK.       =  '  o .~  ^^^  ^      =  ~/~T~ 
m    2A  m    4A 


Since  only  the  shape  defined  by  this  expression  is  desired,  the  factors 
In  front  of  the  summation  sign  v/ere  dropped.   Then,  in  order  to  obtain 
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a  non-dimensional  profile,  each  value  of  deflection  is  divided  by  the 
deflection  at  x  =  A/2,  the  mid-point  of  the  plate.   The  result  is  a 
series  of  shape  coefficients  which  depend  on  B/A,  C/A,  and  x/A. 

The  calculations  for  w/W  for  various  values  of  B/A,  C/A,  and 
x/A  were  made  much  easier  with  the  aid  of  the  digital  computer.   The 
short  program  used  to  make  the  calculations  is  included  at  the  end  of 
this  appendix.   Table  D6  is  a  list  of  the  input  variables  required 
for  the  program.   This  same  program  is  incorporated  into  the  overall 
ice/plate  interaction  program  discussed  in  Appendix  E.   The  output 
from  the  plate  deflection  profile  and  transverse  strain  computer  pro- 
gram is  given  in  Tables  Dl  to  D5.   A  typical  plot  of  w/W  versus  x/A 
is  shown  in  Figure  Dl.   The  variance  in  deflection  profile  for  a  con- 
centrated load  (C/A  =  0.02)  and  a  uniform  load  (C/A  =  1.0)  is  illus- 
trated in  the  figure.   Other  assumed  deflection  profiles  are  also 
plotted  in  order  to  provide  a  comparison  with  other  membrane  solutions. 

The  only  reason  for  assuming  a  deflected  shape  of  the  plate 

is  to  determine  the  average  strain  within  the  plane  of  the  plate.   This 

is  accomplished  by  measuring  the  length  of  the  deflection  profile, 

subtracting  the  original  length  of  plate,  and  dividing  the  difference 

by  the  original  length.   If  L  equals  the  length  of  the  transverse 

deflection  profile,  and  A  is  the  width  of  the  plate,  then  the  average 

transverse  strain  is: 

£   ^   L-A     (in/in)  (D2) 

^     A 

In  order  to  determine  L  from  the  w/W,  x/A  relations,  the 
deflection  profile  was  divided  into  a  series  of  straight  lines.   This 
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©-EXPERIMENTAL    DATA    REF.   El 
B/A  =2.0        q  =  IOOO  pG» 
SB.E  TABLE    Ell 
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PARABOLIC     PROFILE 

CYLINDRICAL   PROFILE! 


FROM  ■)  ^  .^     ,  ^ 

TABLE  D5f  C/A=I.O 
B/A  =2.0  j  C/A  =  0.O2 
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0.4 
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Figure  Dl.   Deflection  Profile  along  the  Transverse  Axis  of  a  Partially 
Loaded  Simply  Supported  Plate 
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linear  approximation  of  the  curve  allowed  rapid  computation  of  L  by 
the  computer.   In  addition,  the  parameter  W/A  must  be  introduced  so 
that  the  strain  calculation  is  related  to  the  actual  plate  deflection. 
Figure  D2  represents  a  certain  segment  of  the  profile  curve.   By  mul- 
tiplying each  segment  of   A  (w/W)  by  W/A,  the  arc  length  A  L  will 
represent  an  actual  segment  of  the  deflection  profile. 


DEFLECTION    PROFILE 
CURVE 


Figure  D2 ,   Element  of  Deflection  Profile  Curve 
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Thus, 


where 


and 

Then,  the  actual  length  of  the  deflection  profile  will  be: 

J 


L  =    /  A  L.  (D4) 

where  J  equals  the  number  of  divisions  of  the  x/A  axis.   With  this 

value  of  L,  the  average  transverse  strain  can  be  found  from  equation  (D2) 

For  the  case  of  a  square  plate,  the  average  strain  in  the  y 
direction  will  equal  the  strain  along  the  x  axis.   If  B/A  >  1.0,  then 
the  two  strains  will  be  different.   The  transverse  strain  will  be  the 
same  as  found  from  equation  (D2).   The  longitudinal  strain  can  be  found 
using  a  similar  relation,  but  the  length  of  the  curve  must  be  adjusted 
to  reflect  the  longer  dimension  of  the  plate.   This  can  be  accomplished 
simply  by  dividing  the  W/A  term  by  B/A  so  that  the  parameter  W/B  is 
used  in  equation  (D3)  in  place  of  W/A.   Then,  if  L .  is  the  length  of 
the  arc  along  the  y-axis,  the  average  longitudinal  strain  is: 

€y  =  Lj  -  B  (D5) 

B 

Finally,  the  principle  strain  at  the  mid-point  of  the  plate 

can  be  determined  from  equation  (ClO)  of  Appendix  C,  which  is  repeated 
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here   for  convenience. 

e-^   ^    €x     "^V'ey  (D6) 

1  -    V  2 

This  is  the  value  of  strain  that  is  presented  in  Tables  Dl  to  D5  and 
Figure  D3  for  various  values  of  B/A,  C/A,  and  W/A.   The  strain  is 
referred  to  as  the  average  transverse  strain  in  the  tables  because 
only  deflections  along  the  x-axis  are  computed  in  the  analysis. 

Also,  it  should  be  noted  that  Poisson's  Ratio  is  equal  to 
0,5  for  the  calculations.   This  is  the  value  that  is  commonly  used 
for  strains  in  the  plastic  range.   Both  1/  =  0.3  and  0.5  were  used  in 
the  calculations  to  test  the  validity  of  the  membrane  formula  and 
1/   =0.5  generally  gave  the  best  results. 
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TRAMSVFRSE  OEFLFCTTOM  PROFILE  AMH.  AVG,  STRMNi 
R/A^  loOOO    ^^Ic,co^,  c  p/^jTn^  0,500 
RATIO  OF  OE^L,  CnEFFICe  AT  X/A  TO  DFFL.o  COEFFIC,  AT  MIP-^T 


X/A  ,50 


C/A 

0  o  0  ?  0 
0,04  0 
0,060 

0,0"  0 
0,1  oo 
0.  ?00 
0,40  0 
0«  ^00 
0,  800 

1  .000 


V!  (  M  A  X  ) 
n, 0003  5 
0. 00141 
Oo  00316 
Oo  00559 
0, 00  36  8 
0,03325 
0, 11500 
Oo  20^91 
0.  2833°^ 
0,3107" 


UOOOO 
1»  OOOO 
U  0000 
1.0000 
1. 0^00 
1«^0  00 
lo  0000 
1. 0000 
1 . 0000 
1.0000 


40 


0  o  o  ]  5  6 
00^162 

0,^171 

Oo  9132 
C. 9  1 Q6 
0o92«l 
Oo  ^^31 
0*9505 
0,9  537 
0. 954^ 


.30 


0,726" 
0.  72  ■5' 8 
Oo  7293 
0,^31^ 
9, 73  3  6 
Oo  "^493 
0  3  7  ^  2  .3 
0.  "05  8 
Oo  31'77 
0  o  3  2  1  3 


:>0 


0»50  34 
0.  '^0^-0 
0,  5  0-'- 9 
0,5  06  1 
0=5077 
0<,5190 
0  ,  5  5  o  p 
0.5307 
0,5^^9 
0  o  6  0  6  3 


10 


Oo  254  6 
0.2  54  9 
0. 2  556 
0,2564 

0.2574 
0o2643 
0, 2"2^ 
0,3021 
0  o  3  1 7  7 
0,3239 


AVFRAGF    TRA'^ISVC'^^F    S  T  R  A  I  r;    FOR        GIVFN    N'AX     Dl^L, 

W/A  0o025  0,050          0,075  0,100          n^i50  0.200 
C/A 

0,020  0,00276  0.,  01  1  02  '  0„  0  ^47^^  0^^4'^6-  0,09668  0c1AS33 

0,040  0,002  76  0,01103    0,^.^^-72  0o0436O  0,00577  n,l6"47 

0,060  0,002  77  0.0110^-    0.02^76  0,04376  0,  0°6  =  0  0,16^'^6^j 

0, 080  0,  00277  0,O|i07    0,0''4c;q  0,043«4  0,0-70^  0,163-<;^ 

0.100  0,00278  0,01110    0,02^-36  0,043q5  0=O':'730  0,1603s 

0,200  0,00283  O.Qii^r    0,02528  O0O4467  0o093-^3  0ol7i36 

0,400  0,002'^5  0,01175    0,  026^^2  0.04647  0,10264  0.17«O8 

0,600  0,00305  0,01?17    ^,02723  0,^A804  0,10591  0,,  13335 

0.800  0,00312  0,012-^5    Oo02""-^5  0.04^^11  0ol':i810  0,18685 

1.000  0,00315  0.01^55     0,07838  0,0^950  0, 10801  0.18813 

Table  Dl 
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TPAN^VFPSe    DFFU:rTICN'    PRTFILF     AMD    AVGa     STRAI 


R/A=     1«?50  PniSSrif-l    S    RATin=    0»S00 

RATIO    OF    OFFL.     CnF^FIC,     AT    X/A    TO     OE^L.,     COFFFIC,     AT     MIO-nT: 


X/A    o^^n 


C/A 

0»0?>0 
0.  OAO 
O.Of  0 
OoO.^O 
0,100 
Oo?CO 
0,400 
0,600 
0«800 
1.000 


W  ( >^  A  Y  ) 
0.  000-'-^? 
Oo  0016Q 
0. 003  3  0 
Oo  006  73 
0,  010^-7 
Oo  040  3  3 

0  o  ?  6  3  2  5 
Oo  3  7  335 
0« 4^112 


loOOCO 
I.  '3  000 
1  <,  0  0  00 
1,00  00 
!♦  O^'OO 
U  00 on 
1. 0^00 
U  00  00 
1,00  00 
U  000  0 


"^'0 


0o921^ 
0  »  0  ?  ^  0 
0  ,  Q  2  3  3 
0. a23  « 
0»  9  2  50 
0  o  0  3  2 1 
Oo  ^446 
Co  ^50  6 
0  o  '^^31 
0.^?3  3 


o30 


0«"'4  06 
0,"^^'14 
Oo  7^'2  7 
0,"7A4''s 
Oo  74  6^^ 
0,73  96 
0o7R76 
Oo  P,06  5 
0  e  -  1  F  6 
0  ,  8 1  7  Q 


20 


0„5  175 
0, 5130 
Oo  5  138 
0  o  5  1  °  - 
0,5212 
0  o  5  :^.  1  0 
0  ,  5  5  "^  8 
0, 5^22 
0  o  5  9  7  0 
0.6012 


10 


0»  2637 
0,2  64  0 
0^  2  64  5 
0  o  2  6  3  2 
0^2661 
o,2  721 
0,2  37  8 
Co  3  036 
^,3156 
0,3  1-^7 


AVf^PAOF    TQ^WSVF'^SF    STRAIN    rOR       01 VM    '-'AX    DFFL 

W/A  0,02  5  0o050  0,0^3  OolOO 

C/A 

0.0^0    0^00246  0„00'^«4  0,0^205  0o0^3^'8  0o0363o  0,15055 

0,040    Oo  00347  Oo  00^35  0o022CA  0,03-^01  ^,0Po''i-6  0,1^067 

0.060    Oo  00^47  0,00OB6  Oo02210  Oo03-'507  n,OR65f^  0,15036 

0o03Q    Oo 00247  0,03^3?  0^^2214  0.03^14  0,08673  0,15111 

0,100    0^00^48  OoOOQOO  Oo0221S  Oo 03^22  0, 03601  0,15141 

0,200    Oo00?52  0,01005  0,02252  0,03^^1  0,0«816  0,13343 

Oo^^OO    Oo  00261  0,010^2  0,0233^  0,^4124  0,0^1  13  0o153^a 

0,600    0,^0^6^  0,0  10""'3  0,02402  0,0'.240  0,  0^362  0,16240 

0.800    0,00274  Oo.0109?  0,02443  0,0^312  0,0^^12  0,16481 

1,000    0,00275  O.OIOQP  0,02456  0,04333  0,09558  0,16555 

Table  D2 
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TRA'lSVEf^SE    OFFLFCTTON    pRFFTLF     /'^■'^'n     A\/Go     STRAIN 
B/A=    1,500  Pni  SSHN     S    PMIT^    0=.  500 

RATIO    OF    OPfL,     r.Q^FFTC.     AT     X/A    TO     DEFl,     COFrFIC,     M'     MIO-PT, 


X/A    ,^0 

»40 

,30 

c  20 

.10 

C/A 

W  (  M  ^  X  ) 

0,0?0 

0„n^0A7 

1 , 00  00 

0.9?A? 

0c'^46R 

0c5239 

0o'67  8 

0,04  0 

Oo  00187 

U  OOOO 

0,c>?4-^ 

00^-^76 

0,5?/^  4 

0,2631 

0.C6  0 

0. C0419 

1»  OPOO 

0of>254 

0,74o.R 

0.5251 

0, 26P6 

0„0R0 

0»  00741 

1,00  00 

0»9?6^ 

Oo  7Cin3 

0,5261 

Oo  2  6*^3 

0»1  0  0 

0. 01153 

1. 00 0  0 

0, ^274 

Oe-^S^l 

P. 5  2  74 

0,2701 

0.200 

0,  <^^^5b 

1, lOon 

0.9339 

0,  76  4'^ 

0, 5  3  6  4 

9,2  756 

0,400 

Oo  1  5^9  6 

1  »  00  00 

0,^453 

0,7900 

r>.,5609 

0  ,  2  Q  P'  0 

0,600 

0,  30376 

1. oono 

0. 9S06 

0  e   '^  0  f  B 

'>  ,  5  3  2  9 

0.3  0''- 2 

O.POO 

0.  43  710 

1,  o-^co 

0,Q5?a 

0,  ^1^-7 

0.5  9-8 

^.^147 

loOOO 

0. 52354 

1,0  0  00 

0o9  534 

Oc  Rl^.6 

0,5991 

0,318  1 

AVFP^G^^  TPxA" 

VJ/^  0.02  3 
C/A 

0,020  0. ^0^30 

0.04  0  0,^0? 3  0 

0,0A^  0, 0023  0 

O.Of^O  0,00731 

0,100  0,002  31 

0,200  0,00235 

0,40  0  '>,002^'3 

0,600  0,00249 

0.^0  0  0,0^2  53 

1,000  0.00254 


■iSv^R^F   STPATN   FOP      oivyr" 

0.05  0         0,0  75  0,100 

O.OC"]"'  0,02'"'55  0,0^634 

0,00^18  0,92  057  0,^'3637 

0.0091^^  0.  0  2  0^^0  0,03  64  2 

0,OOop]  Oo020{3  0,^^64^ 

0.'^0-2^:  0=1^2068  0.03A56 

rs^^^'^o■^f^  0,02098  0,0^70« 

0„ 00^6  8  0,02  169  0,03  83  1 

0,^0 '^94  0,02  22  5  0,03^12  5 

^,0100^  0o0225^  0,03  9'^8 

0,01014  0,0226"  0,04004 


0.150  0,2  00 

0. 9B054  G, 14  035 

0,08  0  60  0,14045 

0, 08  07  0  Oo 14  062 

Oo  0^0^^  0,  14  0^4 

0, 08100  0,14111 

Oc  0^^210  0,1^2  94 

0,0'^  47  3  0,14725 

0,08678  0,15060 

0,^P«0  0  0,15255 

0.08P34  0.15^10 


Table  D3 
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TPANSV^P  ^F     nf=FL-CTTriN    rr^PFTl.r     AND    AVGo     ST'^AP' 


B/A=     lo^f^O  PGISSHM    S    RATin=    0<,500 

RATIO    TF    n^Fl,     C.'llfFK.^     AT    X/A    Tq    DF^L,     CPFFFIC,     AT    ''lO-^T 


X  /  A     .50 


C/A 
0.0?0 
O.OAO 

n.060 

0.080 
0»100 
0,200 
0.400 
0.600 
0«R00 
1*000 


Vi  (MAY) 
0.000^9 
0»  001 96 
Oo  00  44  0 
0.00779 
n. 01?12 
0.  ^^469? 
0,  1  63?0 
0. 32  3^7 
0, ^7052 
0»  5710? 


1 , 00^0 
1  ,  "^0  0  0 
1 . 00  00 
lo 10  00 
1 .  '^  n  n  0 
1,0000 
lo  00^0 
1.  00 no 
1.00 no 
1.00  on 


.^^0 


0,^255 
0,^26  0 
0, 9266 
0.9275 
0  ,  ^"^  ?  -^  5 
n, 0  34^ 
Oo  ^456 
0.  '^  ^0"^ 
0.^^2  7 
0,0  532 


.30 


0  o  7  5  0  ^^ 
0. 75  17 

Oo  7532 

0.7540 

0, 76  65 
C 3 7910 
Oo  906*" 
Oo  R143 
0. ^159 


.20 


0. 5  2  70 
0,5275 

0. 5282 
0.5291 

0.5303 
0 .  5  3  =^  n 

0o5^>24 
9.5332 
0,5952 
0  ,  5  c  >>  7 


.  10 


0.2^-9  8 
Oo2701 
0,2706 
0„ 2  712 
0,  2''2n 
no  2772 
0,2^10 
0,3  04  5 
^o 3143 
0,3173 


AVF'MGF    TRA'^'SV-^SF     STOAJN'    F^!?        GIVFri    MAX     OEFL, 

W/A          0.025  O.O^^n          0on75  OdOO  0,150  0<,  200 

r/A 

0,020  0,00219  0^00876  0,01963  0,03470  0,07638  0,  1^304 

0,0^-<0  O.On.-J-'O  OoOnQ-'T  0,01964  0,  0^-^-73  0,0''6'^4  n«  13^04 

0.060  n.oo?20  0o^^^78  0.01^67  0,n?A77  0,07704  0,13^i29 

0.0^0  O-O-^^PO  0o'^0^7n  0,  91970  0o034"3  0,077i^  0,13440 

0.100  0,00221  0,00881  n, 01974  0,03490  0,07731  n, 134^5 

0.200  0„n0224  O.Onoa^  0,02002  n„n3538  9,07833  0,1363''- 

0,400  0,on?3i  O,00'^?3  0,  n2C6''  0,03652  0,080^4  0,1-^030 

0,600  0.00237  0o0'^a46  0.07  118  0,0^741  0,^8260  0.1^333 

0,800  0,002^1  0,0nc6o  0,  n7]/,P  0-,037^2  0.^536"  0,1^506 

1.000  0.00242  0,00963  0.n7156  no03806  0,0P.396  0,1^551 


Table  D4 
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TRANSVERSE    OEFl-CTinK    prtFILE    AN'C    AVG,     ST'^ATnj 
R/A=     ?,n")0  POISSON    S    «ATin=    Oo  5-0 

RATTH    HF    OZ^L.,     CTF^Fir^     \j    x/A    TH    orrL,     COF-Fir,,     AT    ^MO-nT, 


X/A    o  50 


AO 


0 


?0 


10 


•  C/A 

'^  (  *''  A  X  ) 

0,0?0 

0, 00050 

1,0000 

0,o?6? 

Oe  7513 

0, 5285 

0, 2  70  8 

0.0^^0 

0. 00?01 

1,0000 

0,9266 

0c^520 

Oo  5  2  90 

0,2-M  1 

0.06  0 

Oo  00^h51 

].ooco 

0, 0  273 

0. 7531 

0.5  297 

0,2  715 

0  o  0  .q  0 

0, OOROO 

1  s  0  o  0  0 

O.^PBl 

Oo  7  5A6 

0, 5306 

0,2  72  2 

0,  100 

0, 01 2^4 

U  00  00 

0,'^2<^1 

0,-^563 

0«  531f^ 

0,272^^ 

0..200 

Oo  04B1P, 

1,0000 

0,935? 

0„  7676 

0, ^402 

0,2  781 

OoAOO 

Oo  17315 

1, 0000 

0,045^ 

^•o70i6 

0,5  6^1 

0,2^16 

0.6  00 

0«  33454 

1  o  0  0  0  0 

0,  '^507 

0,  R0"70 

0.5333 

0,^04  7 

0 .  n  0  0 

.'^.  48^6-^ 

1,0  0  00 

0,0  52  6 

0,8141 

0,  59^f9 

0,3141 

uooo 

0,59600 

1 , 00  00 

0.^^31 

0.8156 

0.5977 

0,3170 

AVERAGE    TRANSVERSE    STRAIN    FOR       GIV^N    MAX     OE~L, 
W/A  Oo025  O0O50         0,0''5  OolCO         0.150  O.?00 

C/A 

0,020    0,00213    0,0O84R    OoOT-Ol  0,03361  0,  07^45  '^,]2'^66 

Qo040    0.00213    0,^'^^''-°    0, 01^03  0.,0^>363  0^07450  0,  12^75 

0,060    0,no?l^    Oo  00^=^0    0,01^05  0,03^6^^.  0,0"'4^-  0,12^  =  0 

0,080    0.00213     OcOO^'52     0,01^08  0,03373  0,0^^71  Ool'^OlO 

0,100    0,00 '•'14    Oo  00^5  3    0^0 1-^2  Oo  033  80  0o07^f>5  0,13033 

0,200    0,002]  7    0,00365     0,0193  3  0,03^-26  0,07583  0,1^1^4 

0,400    0,00224    O.OOog-^     Oo  02001  0,035^4  0.07811  0,1356^ 

0.600    0.0022^    Oo  0091^1    0,0204^  0,03617  0o07Q86  0,138E4 

0,800    0,00233     0o039'>p     0»  02077  0,0^,665  0,  08036  0,14015 

1,000    0.00^33    0.0'^93]     0,02084  0,03678  0,0P]12  0,1405o 

Table  D5 
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Table  D6 
LIST  OF  INPUT  VARIABLES  FOR  COMPUTER  PROGRAM  FOR.  PLATE  DEFLECTION 

PROFILE  AND  TRANSVERSE  STRAIN. 

Computer  Definition  or  Symbol  in  Text 

Input  Variables 

BAl  First  value  of  B/A  to  be  used  in  Eqn.  (Dl) 

C2  Increment  of  B/A  to  be  added  to  previous 

value  for  next  set  of  calculations 

C3  Last  value  of  B/A  desired  to  be  used  in 

calculations 

PR  ■•      Poisson's  Ratio 

NI'JA  Number  of  W/A  values  to  be  used  in  calculations 

NXA  Number  of  x/A  values  to  be  displayed  in  output 

WA  W/A 

Output  Variables 

W(l^lAX)  Solution  of  equation  (Dl)  for  x/A  =  0.5, 

the  mid-point  of  the  plate 

x/A  For  plate  dimensions  of  Figure  2 

C/A  For  load  definition  of  Figure  1 

W/A  For  deflection  definition  of  Figure  3 


See  statements  DEFL0004  and  DEFL0006  in  the  program  listing. 
See  Tables  Dl  -  D5. 
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Appendix  E 
COMPUTER  PROGRAM  FOR  ICE /PLATE  INTERACTION 

The  FORTRAN  IV  computer  program  to  calculate  the  mid-point 
deflection  of  a  plate  and  the  amount  of  energy  absorbed  at  the  expected 
point  of  failure  is  listed  at  the  end  of  this  appendix.   The  logic  of 
the  program  follows  the  step-by-step  procedure  given  on  pages  37  to  39. 
The  values  obtained  for  (C/A)   from  step  1  is  used  as  input  to  the 
program.   The  applied  force  F„  and  its  corresponding  W  are  determined 
in  the  first  part  of  the  program.   Then  the  absorbed  energy  U  is  cal- 
culated from  the  strain  relation  at  the  mid-point  of  the  long  side  of 
the  plate.   Next,  the  applied  force  is  increased  by  a  A  F  and  the  pro- 
cess is  repeated  to  find  the  new  value  of  W  and  U.   The  output  format 
of  the  program  is  illustrated  in  Tables  El  to  E19.   The  description 
of  the  input  and  output  variables  and  their  relation  to  the  parameters 
defined  on  the  page  of  Nomenclature ' are  given  in  Table  E20, 

Similar  to  a  hand  calculation  of  W  and  U,  a  constitutive 
relation  between  the  stress  and  the  strain  for  the  plate  material  is 
required  for  the  computer  program.   A  generalized  stress-strain  curve 
was  used  in  the  program  to  obtain  a  stress  for  a  given  total  strain. 
The  general  form  of  the  equation  is: 

or  =  D(  X  +    ep)""    «  D  (e)""  (ei) 

where  6 ^  is  the  plastic  strain  and  G  is  the  total  strain  of  the 
material.   With  the  proper  choice  of  coefficients,  this  equation  pro- 
vides a  reasonable  fit  to  an  actual  tensile  test  stress-strain  curve. 
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It  was  .found  that  by  setting  X  equal  to  the  yield  strain  of  a  given 

material,  and  using  the  bouiidary  condition  (7  =  0     when  G  =  G  ,  the 

u  u 

exponent  n  and  strength  Coefficient  D  can  be  determined.   For  materials 
such  as  steels  which  show  near  constant  stress  for  the  smaller  values 
of  6  ,  the  value  used  for  X  should  be  the  strain  at  the  point  where 
Cf    ceases  to  be  constant.   For  example,  the  stress-strain  curve  for 
mild  steel  exhibits  a  near  constant  stress  at  the  yield  point  up  to 
approximately  6  =  0.02  in/in.   Thus  X  =  0.02  should  be  used  in 
Equation  (El). 

Since  the  generalized  expression  is  always  increasing  with 
increasing  G     ,  the  choice  of  G   is  important  to  ensure  a  proper 
curve  fit.   It  is  usually  difficult  to  determine  from  an  actual  stress- 
strain  curve,  the  actual  value  of  strain  that  corresponds  to   0"  . 
The  best  selection  of  C       is  at  the  point  where  the  (T  -  G   curve  first 
starts  to  level  off.   The  slope  of  the  curve  may  not  actually  reach 
zero  until  just  prior  to  the  failure  strain.   But  the  ultimate  stress 
is  usually  approached  to  within  a  very  small  error  at  a  strain  con- 
siderably smaller  than  the  actual   6  .   Therefore,  in  order  to  main- 
•^  u  ' 

tain  a  good  fit  of  the  curve  betv.^een  G  =  X  and  6  =   6  ,  a  value. 

of   G  u  should  be  used  that  may  be  less  than  the  actual  ultimate  strain. 

From  the  chosen  values  for  X  and   6   ,  the  values  for  D 

and  n  are  determined  from: 

ln(  Q-q/  g^)  ^g2) 

ln(  A/  e^) 


n  = 


and 

Co 


D (E3) 
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The  use  of  the  generalized  stress-strain  curve  allox\fcd  the 
incorporation  of  the  average  transverse  strain  program  discussed  in 
Appendix  D  into  the  overall  program  as  a  subroutine.   The  name  of  this 
subroutine  is  STRAIN  and  its  listing  follows  that  of  the  main  program 
at  the  end  of  this  appendix.   For  an  input  of  the  variables  B/A,  C/A, 
and  W/A  to  STRAIN,  a  value  of  6  j^  as  given  by  equation  (ClO)  is  cal- 
culated and  returned  to  the  main  program.   Then  the  stress  for  6  ■, 
is  found  from  the  generalized  stress-strain  function. 
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F 
(1  RS) 

VI 

{  IN) 

;>05027,  r 

1 , ^0  6  9 

A100?7.,  1 

3. Bl?3 

615077.,  1 

5.2839 

8?00?^, 1 

6. 66? ^ 

0?5C?''.  1 

7-,  9763 

ICE/PL  ATE  INTFRAfTIGN 

ICE  CRUSfilNH  PRESSURH=    45  >?  PS  I 

PLATE  DIMENSIONS  (INCHES) 
TMICKNESS=0 .125 
WIOTH=    54>0  0 
LENGTH^   8  4,00 

MATERIAL  PROPERTIES 

YIELD  ST!lcss=:  56000=0  PST 
ULT>  STRESS=  71200,0  PS  I 
YIELD  STR/\  IN  =  0,0200  IN /IN 
ULT.  STRAINED, 1400  IN/IN 
POISSON  S  RATin=0,500 

AT  f:iD-PT  OF  LONG  SIDE 
S       C/A    STRAIN   STRESS        U 
(LR/IN)  (IN/IN)     (PSI)      (IN-LB) 

7000,0  1,000  0,0230  56967,5  1309o0 

7000,0.  1,000  0,^'^t60  62054.,  5  2851,6 

7577,6  1.000  0.0661  64903.2  4290,  ^^ 

P013.4  1.000  0,0^57  6":'01?.,  ^  5''41.,0 

8366,,  2  loOOO  0,]04B  63701,9  7201,2 


Table  El.   Plate  IB  of  Reference  E3 
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ICE/PL  \]  E     INTERACT  ION 

ICE    CRUSHING    PRESSURE^  43.9    PSI 

PLATE    DI'^'ENSIGMS     (INCHES) 
THICKNESS=0. 182 
WIDTH=  5 A. 0  0 

LENGTH^       8A.00 


F 

(IBS) 


^'.ATERIAL     PROPERTIIS 

YIELD  STRESS-  A1800.0  PSI 
ULT.     STRESS  65500.0    PSI 

YIELD  STRAIK'. 0.0175  IN/IN 
ULT.  STRAIH^O.  1600  IN/IN 
PCISSON    S    RATIO=0.500 


W 
(IN) 


AT  MID-PT  OF  LONG  SIDE 


221810.3  1.8981 

4^3810.3  3.7318 

665810.3  4.9£42 

8878  10.3  6.12  31 


S       C/A  STRAIN  STRESS 

(LB/ IN)  (IN/ IN)  (PSI) 

7637.6  1.000  0.0333  47631.6 

7  742.2  1.000  0.06  60  5  47  27,4 

8696.5  1.000  0.0933  5 870?, 7 

9A39.3  I.OOO  0,1194  61725.4 


11C981C.:       7.1896       10049.3        I.OOO       0.1449       64197.4 


U 

(  IN-LB) 

15  86. 3 
36  12.  8 
5474.4 
73  71.9 
9303  ,7 


Table   E2.       Plaie    2A   of   Reference   E3 
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TCE/PLATR  IMTERACTIPN 
ICF  CRUSHING  PRES?U!^E-   ?41..0  PSI 

PLATE   DP'^^'SIn^!S   (irjciiFS) 

THICKf4ESS-0,?44 
.       WTD^H=  13*50 

LFN1GTH=       21,00 


MATERIAL  PRHPFRTIES 

YIFLO  5;TPESS=  5600030  PSI 
ULT.  STRFSS=  70700»0  PSI 
YIELO  STRA I N= 0.0200  IM/IN 
ULT.  STRA  IN=0, 1400  l"/IM 
PClTSSflN    S    RATin  =  0,500 


AT  MTD-PT  OF  LONG  SIDE 


F 
(LBS) 


W 
(  T  ^' ) 


S 


C/A  STRAIf!       STRESS 

(IN/ IN)  (PSI) 


(  L  R  /  I K! ) 

68323.5       0.3?55       1366490  l.,000  0,0306  53933.  2 

136523.,5       0,6505       13664^0  1*000  O.0612  64028,4 

204723-.5       0,9714       13720.1  UOOO  0,0916  67194,8 

272"23,5       1,2257       14496ol  1,000  0^1186  6?3C7.,7 


U 
(  IN-L8) 

18  05.0 

3  91  8.,  7 

615"^.  3 
8  2 1  ^  >  7 


34112^,5       1,4683       15125,1       1^000       0cl45':^       709f^7>8       10204,8 


Table   E3.       Plate   P2    of   Reference   E2 
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TCE/PLATF     INTERACT  I HM 

ICE    C^.IISHPJG    PRESSURE^       1^0 :0    PS  I 

PLATE     DIMENSIONS     (INCHES) 
THIC1<NESS  =  0-1  19 
WIDTH=  13,50 

LEr)GTH=       21,00 

MATERIAL    PRGPERTICS 

YIELD  STRESS=  56000,0  PSI 
ULT,  STPESS=  70700»0  PSI 
YIELO  STRATN=0o0200  IN/IN 
ULT.  STWAIN=0.1400  IN/UJ 
PGISSON    S    RATI  0^0, 500 

AT    ''':in-PT    OF    L^NG    SIDE 
F  W  S  C/A  STRAIN       STRESS  U 

(LPS)  (IN)  (LR/IN)  (IN/IN)  (PSI)  (IN-LB) 

5^865,0  0.526^  656-^^. 0  KOOO  0->0241  5"7273=4  1383>1 

107665,0  1.0332  6783.9  1,000  0,CA78  62161o7  2971, B 

161''f65>0  1.4329  7336o0  1>000  0.06«^8  64935.4  /■.469,3 

215265,0  l.BO^l  •7-^/i6.A  1,000  0,0?93  66990,1  59.?0>3 

269065,0  2,1707  f^069,4  1,000  0.10^4  6B638.4  7506,0 


Table  E4.      Plate   P3   of  Reference  E2 
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TfE/PLATE     INTER  ACT  I CIM 

ICF    C  RUSHING    PRFSSUf^E=       1?0.0    PS  I 

PLATE    niMFNSTOMS     (IN'CHES) 
TflICKnFSS  =  0,105 
K'inTM=  13.50 

LENGTH^-       21->00 


F 
(LPS) 


MATEf^lAL    PPGPFRTI^S 

YIELD  STPrSS=  ??000.0  p$I 
ULT,  STRrSS=  46100o0  PST 
YIELH  STRAlN=0o0200  IN/IN 
ULT,  STRA  IN=0.,  1600  IN/IN 
PaiSSON'    S    RAT  I  0^0,500 


w 

(  U!) 


340?0  ,  0  0-,  6^92 

680?0,0  1,2120 

102020«0  1,6380 

1360  2  0.0  2 o 030  7 

170020,0  2,4000 


S 


AT  MIO-f^T  OF  LOMG  SIDE 


C/A    STRAIN   STRESS 
(IN/ IN)     (PS  I) 


(LR/IN) 

3360.0  1.000  0.0267  33662.7 

3633->6  1,000  0,0^11  37724.1 

^05^-. 8  1,000  0,0727  40136,2 

4360.6  KOOO  0^0^36  41Q57,1 

4611,8  1>000  0,1140  A3434,5 


U 

( IN-LR) 

8  9^,4 
19?6o3 
2917. 2 
3926.^3 
4950, 2 


Table  E5.      Plate   P4   of  Reference  E2 
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irE/PLATR  INTPP.ACTTHN 

ICE  CRUSHinr;  PRFSSU9E=    87. O  PSI 

PlMf    OIMFNSIGNS  (INCHES) 
THICKNESS=C>06a 
WIDTH=    1, 3.>50 
l.ENGTH=   21,00 

MATERIAL  PROPERTIES 

YIELD  STPf-SS^  32000.0  PSI 
ULT,  STRESS^  41R00>0  PSI 
YIELD  STq/\IM  =  0o0200  IM/IN 
ULT,  STRA  IN=0,1600  IN/IN 
pniSSON  S  RATTn=0=500 

AT  MID-PT  HE  LHMG  SIDE 
F         V/  S       C/A    STRAI^!   STRESS        U 

(IPS)       (IN)     (LB/IN)  ,        (TN/I^1)    (PSI)     (IN-LB) 

24664.-5  0.7379  2176,0  1.000  0^,0193  31B64,3  616,6 

4928^.5  1,3536  2370.4  1,000  0,0370  34627,5  1280,0 

73904.5  1,8709  2571,6  1.000  0,0532  362B5.5  193^^,3 

9852^-,5  2.3545  2724,2  UOOO  0,06-0  37515,4  2536,8 

123144,5  2,8160  2846,9  1,000  0,  0P44  38502,4  324c>,6 


Table  E6.   Plate  P6  of  Reference  E2 
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IfE/PLATE     INTERACTinN 

ICF    CRUSHING    PRFSSURE=       194:.?    PSI 

PLATE    Dl'-'-NSlONS     (INCHES) 
THICKNESS=0.119 
VJT0TH=  13.50 

LEMGTH=       21,00 


F 
(LPS) 


MATERIAL  PROPERTIES 

YIRLD  STRFSS=  42000..- 0  PSI 
ULT,  STRESS^  61000»0  PSI 
YIELD  STRAnvl=0,0200  IM/IM 
lILTo  STP  AIN=0-,1600  IN/IN 
POISSOM  S  RATin-=0>'^00 


W 

( I  r  i ) 


5?2?0.7  0,6R02 

1 044:5  o,''  u  2  40-1 

1566?0>7  1.67A7 

208P20.7  2.0726 

26102  0,7  2,4  442 


AT  MIO-PT  QP  LHNG  SIDE 


S 


C/A  STRAIN       STRESS 

(IN/IN)         (psn 


(Lf\/IN) 

4<?98o0  1.000  0.0312  454^2, R 

5478.4  1,000  0,0505  51074,1 

60RS,3  1,000  0=0846  54410.7 

655q,l  1,000  0,1 OSo  56928.9 

6'-i52,4  1,000  0^1325  58969.9 


U 
(  IN-LB) 

1420,0 

3037o9 

4604:.  6 

6  1  9  S .-,  o 

7813,6 


Table  E7.      Plate   P8   of  Reference  E2 
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IC[:/PI..ATE  INTERACTTHM 

ICF  CUSHIMG  PPESSUPF=   100  0  ^S>] 

PLATF  DTMFNSinNS  (INCHFS) 
THIrKMFSS=0, 188 
WinTH=    13o00 
LFNnTH=   3  6.00 


F 
(IPS) 


M.MERTAL  PRGPFRTIES 

YIELD  STRESS-  38500,0  PST 
UI..T0  STRESS-  59000,0  PSI 
YIELD  STRAIM=^, 0?no  IN/IN 
ULT.  STRA  IM  =  0, 1600  IN/IN 
PGI SSON  S  RAT  10  =  0, 500 


W 
(IN) 


46800,0  0o?843 

140^00.0  0oB5'?4 

23380r^oO  K?766 

3? 73 00,0  1,6?2^ 

4?.0"^00o0  l„c-'44  0 


AT  MID-PT  OF  LONG  SIDE 


S 


C/A    STPAIM   STRESS 
(IM/TM)     (PSI) 


(  L  n  /  I  M  ) 

7?38,0  1„000  0.0?05  38691.0 

7?3P.O  KOOO  0»061A  48472,2 

8053o7  1,000  0„0"6«  5321 60 3 

086t),0  1,000  0,1292  56462,5 

9519,0  1,000  0,1605  5^041,5 


U 
( IN-LB) 

792,7 

297"^,,  2 

5 1  5 1  o  0 

7292,.5 

9479,  1 


Table   E8.      Plate   no.    2    of   Reference   El 
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ICE/PL  Air  I  NT  f'R  ACT  TON 

ICF  C.PU^HING  DRFSSU^r=   100. 0  PS  I 

PLATF     DIHFNSIOMS     (INCIIFS) 
TH1CKMFSS=0. IS^ 
WIDTH=  13.00 

LFNGTH=       }9,?J^ 


MATERIAL    pROderTIES 

YIELD  STRESS^  ^B500oO  PST 
ULTo  STRESS=  nOCOO.O  PST 
YIELD  STPATN=0.0?00  IN/IM 
ULT,  STRAIN1  =  0»1^00  T^!/IM 
POTS  SON    S    RATIO^-O.^jOO 


F 
(LBS) 

?5023,  0 

58  8?  5,0 

976? 5  ;0 

160?? 5.0 
1 ^An?5, 0 
??78?S,0 

?616?5,0 
?n5^!?5-,0 


W 
(  IN) 

0.,  ?^65 
Oo  555R 
0,5^75  2 
lo 105  9 
1 , 30RO 
l.>  5001 
1,68?^ 
U  ^56^ 
?.,  0?5  6 


AT  MIO-PT  OF  LONG  M dE 


S  C/A  STRAIN  STRESS 

(LR/IN)  (IN/IN)     (PSI) 

7?3"..0  UOOO  0-017P  37590. A 

7238.0  1,000  0,0^^13  447^9,8 

72  3  8.0  1,000  0.C6  5C  49175.3 

781  P,  5  l.>0^0  0,086  4  51983,  4 

3  371,6  loOnO  0,1057  54  ir^-,,^, 

3845,6  1,000  0,1 ?46  5604^.7 

9?62,3  1,000  0,1431  57660,8 

9635,4  l,00n  0,1613  59096.? 


t) 

nn-LB) 

6^Q,? 

1^^4.6 
3?M).,  1 
44f'.)r,-8 
57:'  ^,  A 
69ft  1  ,3 

8  2'','),  5 

9  5  \0.  Q 


9974,1       1,000       0.1792       60391.2       lOR.'/.  ,3 


Table  K9 .      Plate   no,    3   of  Refei,-uce  El 
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[CE/PLATF  I^!TER'^CTI^jN 
ICc  CRLSHIN';  PKESSURE-   100.  D  PS  I 

plate:  r: 1 1  uu s i :j ns    (inches) 

THICKNESS^:.  133 
WIOTH^  13.03 

L  E  K  G  r  H=       ]  3  .  '1 0 


F 
(L  BS) 


MATERIAL     PKHPEPvT  I  ES 

YIELD  STRESS^  38  5  0:'.''  PS  I 
ULT.  STRESS^  59000.0  PSI 
YIELO  STf'Ari^O.OZOO  IN/I?^ 
ULT.  ST^,AIi\=0.1600  I,N /!,'■: 
pcjI'SSOm    s  rat  in  =  0.  50  0 


{  IM  ) 


16^- CC:^  C.  17  17 

50700.0  0.5151 

5^5;:'C.C  0.8^82 

UviSOO.O  1.C76A 

1-52  10C.O  1.2  0o5 

1  a  5 0 0 C,  0  1 .  '« ^  {  6 

2  19  7:iO.C  ].'.7  50 


0  645 .2 

2  53  50C.C       1.8  555       10  346.?. 
28730C.O       2.0296       10409.0 


AT    MID-PT    r^f    L^NG    SIOE 

S                 C/\  STRAIN  STRESS                    U 

(  LE./IN  )  {  r^/IN)           C'S]  )  (I  -^-LP 

7233.0  1.000  0.0142  35870.3  508.3 
72  3  0.0  1.000  0.C42  5  4  4  944.9  1910.5 
7325.9  1.000  0.0704  4934^.0  3503.9 
80-31.4  1.000  0.0935  52055.1  4949..2 
3.0  30.0       1.000  0,116  2  5  5247.1  6417.9 

9105.1  1.0^0  0.13P.1  57  24  1.9  7903.8 
l.?:0  n.15^'6  58966.2  9403.3 


l.OnO       0.1S07       60491,2       i:930,C 
l.OOC       0.2015       61062.6        12467.5 


Table   ElO,       Plate   no.    5    of   Reference   El 
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.   ICC/PLATO  INTERACTION 

ICE  CRUSHTMG  PRESSU^,E^   100>0  PSI 

PLATE     DIMENSIONS     (INCHES) 
THICKNESS=0^?50 
WIOTH=  13.00 

LENGTH^       2  6,00 


MATERIAL  PROPERTIES 

YIELD  STRESS^  66000,0  PSI 
ULT,  STRESS=  83000.0  PSI 
YIELD  STRA I N= 0,02  0  0  U!/IN 
ULT.   STRAIN=0ol200  IN/ IN 

p^Isso^sl  s  rati 0=0.500 


AT  MID-PT  OE  LONG  SIDE 


r 
(LPS) 


w 

(  IN) 


C/A  STRAIN       STRESS 

(IN/ IN)  (PSI) 


(LE/IN) 

33R00,0       0,1167       16500,0  1,000  0,0114  61412,0 

101300,0       0o3497       16'500,0  KOOO  0^0341  70669.7 

168P00>0       0,5827       16500,0  1.000  0,0569  75439.6 

236300.0       0.8157       16500,0  1,000  0,0706  78756,4 

303800,0       1,0289       16917,1  1,000  0,1013  81223,2 


U 
( IN-LB) 

69^.4 

2412,0 

42^0.,  5 
6:>7  0,2 
8  2  2^,8 


371300*0       1,2075       17513,8       1.000       0.1212       83107.2       10074.1 


Table  Ell.      Plate   no.    4   of  Reference   El 
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ICr/^LATr     INTERACTION 
ICE    CRUSHING    PRnSSURE=       600,0    RSI 

PLATE  ni MENS  inns   (inches) 

THICKMESS=0.750 
WTDTH=  2 A3  00 

LENGTH^       48^00 


F 
(IPS) 


MATERIAL  PROPERTIES 

YIELD  STRESS^  3P5nOoO  PS  I 
UL.T,  STRFSS=  5^000,0  PSI 
YIELD  STPA  IN  =  0^0200  IN/TN 
ULT,  STRA IN=0, 1600  IN/IN 
POISSON  S  RATTO^0«500 


W 
(IN) 


s 


AT  NID-PT  OF  LONG  SIDE 
C/A    STRAIN   STRESS        U 

(PSI)  (IN-LB) 


(LB/ IN)  (IN/IN) 

423?6o  0'  0.8001  23875,0  0,350  0,0"^5"^  50595, "^  3829,4 

142335,9  1,>2013  23875.0  0,642  0ol3R8  57301.5  7952.3 

242335,^  U3'-^50  28^75,0  0,837  0,1^11  60518-,4  10^583  8 

342335,9  U3636  2B875.0  0»  995  0,215?  62702,8  13495.2 

442335,0  1,5'''2A  2^875,0  1,000  0.2-46  64374,1  15749.1 

542335,9  U7079  28875,0  1,000  0.2709  65735,0  17807,3 

642335,9  1,8308  29316,0  1,000  0,2^?^^  66775,7  19527,6 


Table   E12.      Arbitrary   Plate   no.    1 
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TCE/PLATE     INTERACTIPN 

ir.F    CRUSMIMG    PpfSSURF=       600,0    PS  I 

PLATF    DIMEK'SIONS     (INCHFS) 
TtnCKNESS-^1.000 
WIDTH--^  ?4,00 

LENGTH=       48,00 


F 
(LBS) 


.MATERIAL  PRHPERTIFS 

YIELD  STPESS=  38500,0  f'SI 
ULT.  STRESS=  59000,0  PS  I 
YIELD  STPA IN=0,C?00  IN/IN 
ULT.  STRAIN-0,I  600  If'/ IN 
PGISSHM  S  RAT  I  n:=0.  500 


W 
(IN) 


AT  MID-PT  OF  LOMG  SIDE 
C/A    STRAIN   STRESS        U 

(PS  I)  (I\'-LB) 


(  L  R  /  I M  )  (  I  ^i  /  I  N  ) 

10A5^i3-.9  0,B?59  33500.0  0,550  0,1189  55515.  0  660?,  8 

?04543,>^'  0,9"^56  38500,0  0,769  0,1664  5^4"^4,3  9894,4 

304543,°  UO?!"!  38500,0  0,939  Oo  2030  61954,5  12  576,6 

40454?, 9  U1063  3^500,0  1,000  0,2340  63786,"'  14923,-8 

50'45A3,9  1,2355  38500.0  1,000  0.2613  65249,5  1704Ro8 

6045^13,^  1.3524  33500,0  l.,000  0,2860  66471,7  19011,6 


Table  E13.   Arbitrary  Plate  no.  2 
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ICP/PLATf:     INTFPACTIOM 

ICE    CfUJ.SHING    PRFSSURF=       600.0    RSI 

PLATH     DIMFr^SIGNS    (INCHES) 
THICKNESS=1>250 
WIOTH=         24,00 
LF.MGTH=       4  8a  00 


F 
(LPS) 


MATERIAL  PROPERTIES 

YIELD  STRESS=  38500.0  PSI 
ULT,  STRFSS=  59000,0  PSI 
YIELD  STRA  IM  =  0.0200  IN/ir.' 
ULT,  STRAIN=0, 1600  IN/ IN 
POTSSON  S  RATn  =  0^500 


W 
(  IN) 


S 


AT  iMTD-PT  or  LONG  SIDE 
C/A    STRAIN   STRESS        U 

(IN/IN)     (PSn     (IN-LR) 


(LB/ IN) 

279^35.7       0,812^      48125„0  0.,900  0.1946  61421.0  nc)54.,0 

^7993  5.7       0>R577       43125»0  1,000  0,2267  63377.>  1  14369,9 

479935.7       0.9640       48125.0  1..000  0.254f^  64915^5  16542,7 

579035,7       L,  05O-T       AR125,0  1,000  0.  2P^01  6618R,R  185^1^4 

6799^5,7       1,]'^74       48125,0  1,000  0^3033  6727B..  5  20406,9 


Table   E14.      Arbitrary   Plate   no.    3 
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ICF/PLATE     INTERACTION- 

ICF    CRUSHING    PRFSSURF.=       600.0    PS  I 

PLATE    DIMFNSIONS     (INCHES) 
THICKNESS-1.500 
VlinTH=  24«00 

LENGTH-       4  8,00 


MATERIAL    PROPERTIES 

YTELn  STRESS-  3R500>0  PSI 
ULT,  STPFSS-  59000.0  PSI 
YIELD  STRAIN=0a0?00  IN/IN 
ULT,  STRAIN  =  0<.1600  IN/IM 
POISSON    S    RATin-O. 500 


F 

( L  n  s ) 


w  s 

(IN)  (LB/  IN) 


AT  f'in-PT  HE  LONG  SIDE 
C/A    STRAIN   STRESS        U 

(IN/IN)     (PSI)     (IN-LB) 


6oi?00o0  0.6817  57750^0  1,000  0.2162  62764,0  1-^5^2.6 

701200.0  U0314  57750.,  0  1,000  0,3272  63333,2  22358,5 

P9120O0O  1.>0QA7  57750.0  UOOO  0,3473  6^173..  1  24021,1 

991200,0  1,1545  5775O0O  1,000  0.,  3662  69932,4  25611,0 


Table  E15.   Arbitrary  Plate  no.  4 
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TCF/PLATl^     IMTFRAC.TinN 

ICE    CRUSHING    PRFSSURF=       600.0    PS  I 

PLATF     Dlf'FNSIONJS     (INCHES) 
THICKMPSS=0.750 
V/IDTH=  24,00 

LENGTH---       4 80  00 


F 
(LPS) 


MATERIAL  PPGprRTTES 

YIELD  STRE$S=  66000.0  PSI 
ULT,  STRESS^  83000,0  PST 
YI^LO  STRA IN=0,0200  IN/IN 
ULT,  STRAIM=0,T200  IM/IN 
POISSON  S  RATI0=0o500 


w 

(  IN) 


AT  MIO-PT  HF  LHNG  SIDE 
C/A    STRAIN   STRESS        U 

(IN/TN)     (PSI)      (IN-LB) 


(  L  B  /  I N  ) 

104543,9  0.6423  49500,0  0,550  0.06^4  77382,3  5368.8 

254543.'^^  0,'^a33  49500>0  0,35v8  0,1003  8]"14.0  8^67,9 

404543,9  0,8605  4^500,0  1,000  0ol365  84377,4  11515,8 

554543,9  1,0074  4^500,0  1,000  0,1598  860^6,6  13"757>4 

704543>9  l.,1^55  4950^.0  1,000  0,1801  8^425,0  15746,1 

854543>9  1.2506  49500.,0  1,000  0,1984  88510,^  17556>9 


Table   E16.      Arbitrary   Plate   no.    5 
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ICE/PLATE    INTERACTION 

ICE    CRUSHING    PRESSURE^      600.0     PS  I 

PLATE    CIMENSIONS    (INCHES) 
THICKNESS=I .000 
WIDTH=  24.00 

LENGTH=       48.00 


F 
(LBS) 


MATERIAL    PROPERTIES 

YIELD  STRESS=  66000.0  PSI 
ULT.  STRESS=  83000.0  PSI 
YIELD  STRAIN=0.0200  IN/IN 
ULT,  STRAIN=0.1200  IN/IN 
POISSQN    S    RATIO=0.500 


W 
(IN) 


AT    MID-PT    OF    LONG    SIDE 
C /A  STRAIN       STRESS  U 

(IN/IN)  (PSI  )  (  IN-LB) 


(L8/ IN) 

221183.9  0.5765  66000.0  0.800  0.1009  8I18I.3  8192.5 

371183.9  0.6182  660CO.O  I. 000  0.1307  83914.2  1CS7C.2 

521183. 9  0.7325  66000.0  1.000  0.1549  85755.7  13284.4 

671183.9  0.8313  66000.0  1.000  0,1758  87154.2  15321.2 

.82]183o9  0.9195  66000.0  1.000  0.1944  88285.8  17167. C 

971183.9  0.9999  66000.0  1,000  0.2115  89238.2  18370.6 


Table   E17,      Arbitrary   Plate    no.    6 
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ICE/P1.ATF     I^JT^RACTln^J 
ICE    CRUSHING    PRESS'JPff=       600,0    PS  I 

PLATR    nTMFNSIOMS     (PvlCHES) 
TtnCKNESS^l  ,?50 
WIOTH=  2  4,00 

LENGTH^       4Ro00 


F 
(LPS) 


MATERIAL  pRODFRTIES 

YIFLO  STRPSS=  66000,0  PSI 
ULT„  STRFSS=  83000,0  PSI 
YIELD  STRATM=0,0?00  IN/IN 
ULT,  STRAIN=0. 1?00  IN/ IN 
PGISSOM  S  P AT  I  0  =  0. >  300 


W 
(IN) 


S 


AT  MID-PT  OF  IHNG  SIOE 
C/A    STPAIN   STRESS        U 

(  IN/H')     (PSI  )     {  IN-LB) 


(LP-/ IN) 

6^1?0n..0       O,^^??       R?500,0  loOOO  0.1261  e3531o9  10537,1 

R91?00>0       Oa 7663       82500,0  1,000  0,2^26  8B7AQ,0  17977,7 

lO'"-^12no,0       0,847^       82500,0  loOOO  0,2241  89905,"'  20152,2 

12^1200,0       0,^224       82500,0  1,000  0,2438  ^0878, 6  22158.6 

1491200.0       1,0205       82500..0  1,000  0.2721  92164,9  25082,4 


Table   E18.      Arbitrary   Plate   no.    7 


116 


ICE/PLATF  I MTF PACTION 

ICE  CRUSHING  PRESSURE^   600.0  PST 

PLATE     ^IMENSI■^^!S     (INCHES) 
THICKNESS=],500 
Viir)TH=  2 '1,0  0 

LENGTH^       4  8,00 


F 

(LPS) 


MATE^ML     P'-^QPcRTTES 

YIELO  S^RESS=  66000.0  f>SI 
ULTo  STRESS^  33000,0  ^SI 
YIELO  STRAIN=0o0?00  IN/IN 
UL^o  STRA IN=0,1200  IN/IN 
POTSSCN    S    RATin=0.500 


W 
(IN) 


S 


AT  MID-PT  OF  LONG  SIOE 
C/A    STRAIN   STRESS        U 

(PS  I)     (IN-LR) 


(Lf^/IN)  (IN/IN) 

6^1700.0       0,3977       Q9000.0  UOOn  0>]261  03531, ^^  1053'', 1 

891200.0       0,6386       99000,0  1,000  0,2026  8874^^,0  17Q77,7 

1091200.0       0.7066       90000,0  1,000  0o2241  89905.7  20152,2 

1291200.0       0,7686       9QOOO,0  1,000  0,2^38  90878.6  22158,6 

1491200.0       0.P.579       Q^OOO,  0  1,000  0.2721  92164.0  25082,4 


Table  E19.      Arbitrary  Plate  no. 
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O^o 

U 

x 

in 

equation 

(El) 

6 

u 

Po 

isson's   Ratio 

AF 

Table  E20 
LIST  OF  INPUT  VARIABLES  FOR  THE  COMPUTER  PROGRAM 
FOR  PLATE  DEFLECTION /LOAD  RELATION 

Computer  Definition  or  Symbol  in  Text 

•k 

Input  Variables^ 
YIELD 
ULT 
EY 
EF 
PR 
DF 

FLIM  Last  value  of  F  to  be  used  in  the 

calculations 

H  H 

A  A 

B  B 

CA  (C/A)^ 

PC  P^ 

Q  k 

QD  kd 

Output  Variables    (See  Tables  El  to  E19) 

All  of  the  variables  displayed  in  the  output  are  self  explan- 
atory or  of  the  same  notation  as  used  in  the  text.   The  STRAIN  at  the 
mid-point  of  the  long  side  is  equal  to  Q,  ^.    in  equation  (19)  of  the 
text.   The  STRESS  at  the  same  point  is  equal  to  G .. 


See  statements  FAIL0003  and  FAIL0005  in  the  program  listing, 
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